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SUMMARY 

Donnell  type  equilibrium  and  stability  equations  are  derived  for  stiffened  thin  conical  shells. 

The  stiffeners  are  considered  closely  spaced  and  are  therefore  assumed  to  be  ‘‘distributed”  over  the 
whole  surface  of  the  shell.  In  the  proposed  theory  the  stiffeners  and  their  spacing  may  vary  in  any 
prescribed  manner,  but  here  only  equally  spaced  stiffeners  are  dealt  with.  The  force  —  and  moment  — 
strain  relations  of  the  combined  stiffener-sheet  cross-section  are  determined  by  the  assumption  of 
identical  normal  strains  at  the  contact  surface  of  stiffener  and  sheet. 

The  stability  equations  are  solved  for  general  instability  under  hydrostatic  pressure  by  the  method 
of  virtual  displacements.  The  solution  used  earlier  for  unstiffened  conical  shells,  which  satisfies  some 
of  the  boundary  conditions  of  simple  supports  only  approximately,  is  again  applied  here.  The  effect  of 
this  incomplete  compliance  with  boundary  conditions  is  shown  to  be  negligible  by  consideration  of 
“boundary  work” .  The  solution  proposed  for  stiffened  conical  shells  involves  the  concepts  of  “correc¬ 
ting  coefficients”  and  minimization  of  corresponding  “error  loads”. 

Typical  examples  are  analysed  and  the  effect  of  eccentricity  of  stiffeners  is  investigated. 
Simplified  approximate  formulae  for  the  critical  pressure  of  frame-stiffened  conical  shells  are  also 
proposed. 
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distance  of  the  top  of  a  truncated  cone  from  the  vertex,  along  a  generator 
(see  Fig.  1). 

distance  between  the  rings  (see  Fig.  1). 
defined  by  Eqs.  (97). 

real  displacement  coefficients  defined  by  Eq.  (67). 
expressions  defined  by  Eqs.  (96). 

cross-sectional  area  of  stringer  or  frame  (ring  stiffener),  respectively, 
complex  displacement  coefficients  defined  by  Eqs.  (46). 
distance  between  the  stringers  at  x  =  1  (see  Fig.  1). 
npmv/£  ,  see  Eq.  (161). 

real  displacement  coefficient  defined  by  Eqs.  (46). 
defined  by  Eqs.  (76)  to  (78). 

Eh3/12  (1  -  v2) 

distance  of  the  centroid  of  the  stringer  cross  section  from  the  shell 
middle  surface  (see  Fig.  1). 

distance  of  the  centroid  of  the  ring  cross  section  from  the  shell  middle 
surface  (see  Fig.  1). 
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=  moduli  of  elasticity  of  sheet  or  stiffeners,  respectively. 

=  effective  moduli  of  elasticity  of  stiffeners  (see  Section  16). 

=  expressions  defined  by  Eqs.  (128)  to  (142). 

=  defined  by  Eq.  ( 1  IS). 

=  Pct/Pcr»  defined  by  Eq.  (159). 

=  shear  moduli  of  the  stiffeners. 

=  thickness  of  shell. 

=  moment  of  inertia  of  stringer  or  frame  cross-section  respectively,  about 
the  line  of  reference  (the  middle  line  of  the  sheet). 

=  moment  of  inertia  of  stringer  or  frame  cross-section,  respectively,  about 
their  centroidal  axis. 

=  torsion  constants  of  stiffener  cross-section. 

=  expressions  defined  by  Eqs.  (98). 

=  correcting  coefficients  defined  by  Eqs.  (87)  to  (92). 

=  spring  constants  defined  by  Eqs.  (106). 

-  a  (x2  -  1),  see  Eqs.  (161). 

=  integers. 

=  moments  and  forces  prior  to  buckling. 

=  additional  moments  and  forces  caused  by  buckling . 
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=  total  moments  and  forces  during  buckling. 

=  moments  and  forces  acting  at  the  boundaries  prior  to  buckling. 

=  number  of  displacement  terms  . 

=  expressions  defined  by  F.qs.  (69). 

=  hydrostatic  pressure,  critical  pressure  of  stiffened  conical  shell  and 
critical  pressure  of  unstiffened  conical  shell,  respectively. 

=  critical  pressure  of  equivalent  cylindrical  shell , 

=  external  loads  . 

=  radii  of  small  or  large  end  of  truncated  cone,  respectively. 

=  complex  number  ,  s  =  y  +  i  n /3 
2  2 

=  tQ  +  cQ  ,  see  F.q.  (163)  . 

=  number  of  circumferential  waves  . 

=  t/cos  a,  see  Eq.  (161) 

=  expression  defined  by  F.qs.  (117)  and  (126). 

=  displacement  along  a  generator 

=  non-dimensional  displacement  along  a  generator  =  u*/a 
=  total  potential  energy  prior  to  buckling  . 
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=  circumferential  displacement 

=  non-dimensional  circumferential  displacement  =  v*/a 
=  radial  displacement. 

=  non-dimensional  radial  displacement  =  w*/a 
=  axial  coordinate  along  a  generator, 
s*  non-dimensional  axial  coordinate  =  x*/a 

=  ratio  of  the  distance  of  the  bottom  of  a  truncated  cone  from  the  vertex, 
to  that  of  the  top , 

=  radial  coordinate . 

=  distance  of  the  centroid  of  the  stringer-shell,  or  ring-shell  combination 
from  the  middle  surface  (see  Fig.  1). 

=  cone  angle 

=  rr/lnex2,  see  Eq.  (SO). 

=  defined  by  Eq.  (56),  =  (1/2)  [1  -  v/(\  +  *J01)] 

=  middle  surface  strains. 

=  Fi\eia/boD 
=  E2A2e2a/a0D 

=  bi*oi//bo^ 
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=  changes  of  curvature  and  twist  of  the  middle  surface  . 

=  non-dimensional  changes  of  curvature  and  twist  of  the  middle  surface 
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=  Poisson’s  ratio . 
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=  (1- v2)  E2A2e2/EaQha 


Subscripts  following  a  comma  indicate  differentiation. 
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I.  INTRODUCTION 

In  order  to  increase  the  resistance  of  shells  to  buckling,  they  are  strengthened  by  stiffeners.  In 
this  manner,  the  critical  load  can  be  increased  several  times  by  only  little  addition  of  material.  It  is 
assumed  that  the  buckling  is  of  the  general  instability  type,  that  is,  the  shell  and  its  stiffeners  buckle 
together. 

The  stiffeners  are  considered  closely  spaced,  and  are  therefore  assumed  to  be  “distributed”  over 
the  whole  surface  of  the  shell.  In  the  proposed  theory,  the  stiffeners  need  not  be  equal  and  equally 
spaced,  but  may  change  in  any  prescribed  manner.  The  present  report,  however,  deals  with  conical  shells 
stiffened  by  equal  and  equally  spaced  frames  (rings)  and  stringers.  This  is  the  usual  way  of  stiffening, 
though  not  necessarily  the  optimal  one.  There  may  be  some  other  law  of  stiffener  distribution  which 
would  yield  maximum  stiffening  for  a  given  addition  of  material.  This  optimization  problem  is  not  con¬ 
sidered  here,  but  could  be  solved  by  the  proposed  method  of  solution. 

The  relations  between  the  strains  and  the  internal  forces  and  moments  of  the  combined  stiffener- 
sheet  section  are  found  by  the  assumption  that  the  normal  strains,  in  the  stiffener  and  in  the  sheet,  are 
equal  at  their  point  of  contact.  Thus,  the  eccentricity  of  the  stiffeners  relative  to  the  sheet  is  taken 
into  account.  The  analysis  permits  the  sheet  and  the  stiffeners  to  be  of  different  materials. 

The  middle  surface  of  the  shell  (without  stiffeners)  is  taken  as  the  surface  of  reference.  The  stress- 
strain  relations  in  this  surface  are  assumed  as  in  an  unstiffened  shell. 

It  is  assumed  that  for  general  instability  a  stiffened  conical  shell  buckles  in  a  mode  similar  to  that 
of  an  unstiffened  conical  shell.  Hence,  the  displacements  which  were  used  for  unstiffened  conical  shells 
(Ref.  1)  could  be  applied  to  the  problem  of  buckling  of  a  stiffened  conical  shell.  These  displacements 
satisfy  some  of  the  assumed  boundary  conditions  of  simple  supports  only  approximately,  and  imply  fic- 
titions  elastic  restraints.  The  effect  of  these  restraints  was,  however,  shown  to  be  negligible  (Ref.  2). 
Here,  another  method  is  proposed  to  estimate  the  effect  of  the  partial  fulfilment  only  of  these  boundary 
conditions  upon  the  value  of  the  critical  load.  The  method  is  based  on  consideration  of  the  “boundary- 
work"  -  the  work  done  by  the  internal  forces  and  moments  at  the  boundaries. 
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The  method  of  virtual  displacements,  used  for  solution  of  the  present  problem,  stems  from  the 
principle  of  virtual  work.  It  only  requires  that  the  displacements  fulfil  geometrical  boundary  conditions. 
The  displacements  used  here  fulfil  rigorously  the  boundary  conditions  of  zero  radial  displacement,  while 
fulfilling  only  approximately  the  boundary  condition  of  zero  displacement  in  the  circumferential  direction. 
A  method  of  satisfying  this  boundary  condition  in  general,  although  not  for  ifivery  term  of  the  displacement 
series,  is  then  proposed. 

For  unstiffened  conical  shells  the  displacements  solve  the  first  two  stability  equations  exactly. 

In  the  case  of  stiffened  conical  shells,  however,  these  equations  are  not  solved  exactly  by  them.  This 
occurs  on  account  of  the  additional  terms  introduced  by  the  frames  and  stringers.  But  the  same  displace¬ 
ments  may  be  used  in  the  following  manner.  In  the  first  two  stability  equations,  the  terms  which  do  not 
lend  themselves  to  solution,  are  replaced  by  terms  which  do.  The  latter  are  multiplied  by  coefficients 
called  “correcting  coefficients”.  Thus,  “corrected”  stability  equations,  which  can  be  solved  exactly  by 
the  displacements,  are  obtained.  The  original  terms,  removed  from  the  first  stability  equations,  and  the 
“correcting”  terms,  which  replace  them  there,  with  opposite  sign,  are  added.  These  sums  are  called 
“error-loads”.  The  “correcting  coefficients”  are  calculated  by  equating  the  virtual  work  done  by  the 
“error-load”  to  zero. 

In  the  analysis  the  “effective  sheet  length”  is  considered  as  a  reduction  in  the  moduli  of  elasticity 
of  the  stiffeners. 

A  simpler  approximate  method  of  calculation  is  obtained  by  neglecting  the  eccentricity  of  the  stif¬ 
feners.  Then,  “correcting  coefficients"  are  not-needed,  and  the  calculations  become  easier. 

A  simple  approximate  formula  for  calculation  of  the  critic  1  pressure  of  a  frame-stiffened  conical 
shell,  by  consideration  of  an  equivalent  cylindrical  shell,  is  proposed.  This  formula  is  based  on  a  simi¬ 
lar  one  for  unstiffened  conical  shells. 

Some  typical  cases  are  calculated  by  the  above  methods.  It  is  shown  that  the  effect  of  the  “boun¬ 
dary  work”  upon  the  value  of  the  critical  load  is  small.  Frames  (rings)  increase  the  resistance  of  the 
shell  against  buckling,  under  hydrostatic  pressure,  considerably.  It  is  shown  that  the  placing  of  the 
frames  is  of  importance.  Frames  on  the  inside  of  the  shell  yield  higher  general  instability  loads  than 
frames  on  the  outside.  Stringers  are  much  less  effective  in  stiffening  of  shells  under  hydrostatic  pressure, 
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and  the  effect  of  their  eccentricity  is  opposite;  outside  stringers  yield  higher  critical  loads  than  inside 
stringers. 


2.  PREBUCKLING  EQUILIBRIUM 


The  strain-displacement  relations  at  the  middle  surface  of  a  deformed  conical  shell  used  in  the 
derivation  are  those  given  by  Love  (Ref.  3). 

f  =  u 


((f>  =  v,<f>/x  s*n  a  +  U//x  “  w  cot  a/ x 

y*</,  =  v,*  -  v/*  +  v/x  8*n  a 

The  curvatures  are  defined  as 

“  w,« 

V  *  w,,/*  +  w,^/x2  ain2  « 

KX<t>  m  W,x^/X  ain  ®  -  W,^/* 2  sia  « 


(1) 


These  are  the  curvature  displacement  relations  obtained  by  Seide  (Ref.  4)  by  omitting  the  terms  involving 
the  circumferential  displacement  v  from  Love’s  definition,  on  account  of  their  negligible  effect  in  cylin¬ 
drical  shells  and  vanishing  in  the  case  of  circular  plates. 

The  analysis  is  written  in  non-dimensional  form,  and  the  non-dimensional  distances,  displacements 
and  curvatures  are  defined  by  the  equations  : 
x  «*  x*/a 


z  i 

u 

V 

w 


z*/a 

u*/b 

v*/a 
■  w*/  a 


** 


a  k_ 


a  k 
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K*<t> =  a 


(3) 
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Virtual  displacements  are  applied  to  a  segment  of  the  shell,  which  is  in  equilibrium.  Hence,  by 
the  principle  of  virtual  work  (the  virtual  work  done  by  the  stresses  must  be  equal  to  the  virtual  work 
done  by  the  external  forces  and  moments), 

0=SUo=/2  /  [N„&s  +  N^0S^  +  N^0Sy^  -M,0S(K,/a)-M^0S(v/a) 

4>\  *1 

+  ^os(^/a)-‘VS(V/a)  "  (q*a5u 

+  q^aSv  +  qx  a  S  w)]  a2  x  sin  a  dx  d  0 

4*2  — 

/  lxa[Ni0  5(au)+  Nx<^0  5(av)-  Mi0  5(W>I)  +  Mx^0  5(wf^)/x  Sina 
4*i 

x~x  2 

+  Qx0  S(aw)]  j  sin  a  d  <f> 

I==xi 

“  J[N^05(av)  +  N^x0S(au)  -M^0S(w^)/xsina-M^x0S(wfS) 


4>=4*  2 

+  Q  s(aw)  ]  ad x 


where  the  index  zero  indicates  the  state  of  the  shell  before  buckling,  and  the  barred  quantities  are  the 
external  forces  and  moments  acting  on  the  boundaries. 

It  is  assumed  that  Nx<^  is  carried  only  by  the  sheet  (i.e.  the  stiffeners  do  not  transmit  sheer). 
Hence,  it  is  assumed  as  in  Ref.  4  that 

The  work  done  by  the  internal  forces  Qx  and  Qy  is  neglected  in  the  expression  of  the  virtual 
work,  as  the  theory  developed  is  a  Donnell  type  theory. 

Substitution  of  Eqs.  (1)  and  (2)  into  Eq.  (4),  and  integration  by  parts  yields ; 

0  -  SU0  «  -  /  f2  1  [(x  Nx0)>x/a  x  -  /ax  +  Nx^/ax  sin  a  +  qx]5u 
4>i  xi 
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02 
+  / 

01 


+  [Ns6o,^/aX8in«  +  (x2Nx94o),,/ax2+ 

+  [(x  M«oU/a2)C-M^o.«/a2x  +M9io.^/a2x2  si“2« 

~^x  ^x</iO  &2  x2  Sln  a  +  (x  M^6x  o\x<j>/ &2  x2  ®‘n  a 

+  N<£o  cot  “  /ax  +  qJSw  |  a3x  ain  a  dxd<£ 

I  ax[(Ni0  -r  N^q)  a  5u  +  (N^0  -  N^0)aSv  -  (Mx0  -  Mx0 )  S  ( w>x )] 
+  [U  Mx0)>x  -  sin  a  +  M^x0  sin  „-ax  Qx0 

+  sin  sin  a  &<f> 

X  X  , 


X/'(N0O  -  N^o>aSv  +(N^0-N^0)aSu-(M^0-M^0)S(w^)/x8ina 


+  fM^o.^/xa8ina-(xMx<?io)>x/ax+  (xM(?Sx0)>x/ax  - 
0=02 

“  Mfo,,/a  ]a5wj  adx 
0=01 

+  %-^0)a«wT2V2  -  <M0*o  ~  Vo)aSw  ,“V2 

X-Xj  0-0  j  X=Xj  0=0j 


Eq.  (6)  yields  the  following  equilibrium  equations  and  boundary  conditions  : 
In  the  shell 


(xNxo),x/ax  ~  N<£0/ax+  Nx0o,0/ax  sln  a  +  qx  =  0 

N0o.0/axsln«  +  (*2  N„  0  o\x/a  x2  +  ^  -  0 

<XMxo),xx/a2x  -  M0O.x/a2x  +  V.00/a2x2  B“*“ 

“  (x  M*0o\x0/a2  x2  ain  «  +  <x  M0,o>,«0/a*  x2  sIn  «  +  N^0  cota/ax +  qx  =0 
Along  the  circles  x  =  Xj  and  x  =  x2 


(6) 


(7) 

(8) 

(9) 
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N*o  -  N*o 


N,0o  “  Nx0.o 


M  n  =  M  n 

xO  xO 


(xM,olx  -  %0  -  M^0f^/sin  a  +  M^x0(^/8ina  =  ax  Qk(J  -  M^^/sin  a 
Along  the  generators  0  =  0X  and  0  «=  02 


N0o  =  N0o 


M 


^0x  0  V  0 


and  at  the  comers  of  the  segment 


^x0O  ^x0o 


M  ,  =  M  , 

0x0  0x0 


(10) 

(H) 


M0O  -  M0O 

(12) 

/*x  =  Q0O  +  Vo.x/a 

(13) 

(14) 


The  geometrical  boundary  conditions  are  not  discussed  here,  as  it  has  been  assumed  in  the  derivation 
that  the  forces  acting  on  the  boundaries  obey  any  given  geometrical  boundary  conditions. 


3.  EQUILIBRIUM  DURING  BUCKLING 

The  equilibrium  at  buckling  is  obtained  by  consideration  of  the  additional  virtual  work  during  buck¬ 
ling.  The  displacements  are  now  the  additional  displacements  caused  by  buckling,  and  the  prebuckling 
displacements  are  assumed  to  be  small  so  that  also  the  additional  displacements  can  be  related  to  the 
undeformed  geometry  of  the  shell. 

The  stretching  of  the  middle  surface  introduces  the  following  nonlinear  terms,  which  have  to  be 
added  to  the  strains  of  Eqs.  (1), 

=  (w,x>2/2 

*0  =  (w,0>2/2  x2  sin2  «  (IS) 

Vx0  “  W,x  W,0/X  8in  « 
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The  connection  between  forces  and  moments  prior  and  during  buckling  is 
N*  =  N0  +  N 
M'  =  M0  +  M 


(16) 


Where  N  and  M  are  the  total  forces  and  moments  during  buckling,  NQ  and  MQ  are  those 
prior  to  buckling,  and  N  and  M  are  the  additional  forces  and  moments  caused  by  buckling. 

Now,  if  to  the  shell  is  given  virtual  displacements,  the  virtual  work  8U  done  during  buckling 
must  also  vanish,  since  the  shell  is  in  a  state  of  equilibrium.  The  internal  forces  and  moments  NQ 
and  in  Eqs.  (4)  must  be  replaced  by  N  and  M  from  Eqs.  (16),  and  the  nonlinear  virtual  work 
done  by  the  membrane  forces  prior  to  buckling  Nx  Q  ,  N^  and  Nx(^0  , 

ASU  =  /  /  (NxQ  S  (’x  +  N^0  +  Nr(^o  8  y’x,)  x  a2  sin  adxd0  (17) 

must  be  added. 

After  Eqs.  (7)  to  (14)  are  also  taken  into  account,  the  expression  for  the  virtual  work  during  buck¬ 
ling  becomes  : 


0 


.2  2 


0  -  SU  =  -  /  /  1  [(xNx  )jX/ax  -  N^/ax  +-  N^^/a  x  sin  a]  Su  +  [  N^^/ax  sin  a, 
0!  *1 

+  (x2N^U/a-2lSv  +  [(x  Mx)>xx/a2x  -  M^x/a2x 
+  %,4>4>/*2  x2  8i"2  a  “  (x  V,*^2  x2  sin  «  +  (x  M^x\*0/a2  x2  sin  « 

+  cot  a/a  x  +  (x  Nx0  w>x)>x/ax  +  (N^0  w^)^/ax2  sin2  a 

+  ^xcfco  sin  a  +  (N*<£o  /ax  sin  a]  Sw  la3  sin  a  dxd<^> 

02 

+  /  1  ax  [Nx  aSu  +  Nx<?i  aSv  -  MxS(w>x)]  +  [(xMi)jj  -  -  M^^/sin  a 

01 


(18) 


+  M^jc,q6 /sin  a  +  Nx0  axw^  4  NX(;6o  aw>(jS/sin  a]  aSw  |  sin  a  d  <f> 
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+  ¥  I  Nj  aSv  +  N^x  aSu  -  M^Mw  6)/x  sin  a  +  [M^  ^/ax  sin  a 
xi 


~  (xMx^,x/ax  +  (x  M<f>x\x/*X  +  w,9S/x  sin  “  +  Nx^0  w,x^a5w  •  2  adx 

x=x2  <t>=<f>2  x=x2  c£=<£2  ^  ^ 

+  aSw  I  I  -  aSw  I  I 

*=*1  ^=01  *=x!  <£=<£i 

Hence,  the  following  stability  equations  and  boundary  conditions  are  obtained: 

(xNx),x/ax  “  N^/ a x  +  *x  sin  a  =  0 

N^/ax  sin  a  +  (x2  Nx0),9i/ax2  =  0 

(xM,),,./a2x  -  ^<£,x  / fl2  x  +  M0,^/a2  *2  3in2  tt-(*M,^)^/a2x2  sin  a 
+  (x  V„^/a2x2  sIn  a  +  N<£  cot  a/ax  +  (x  N*o  W,x),x/aX 

+  °V>  W,^),^/ax2  sin2  «  +  (Nx<£o  W,x),«^/aX  sin  “  +  <Nx<£o  w,0\,/axsin  «  =  0 

Along  x  =  Xj  and  x  =  x2 


Nx  =  0  or  u  =  0 


Nxi^  =0  or  v  =  0 

(x  Mx)>x  -  -  Mx^/sin  a  +  M^/sin  a  +  NxQ  xa  w>x  +  aNx^  w^/ein  a  =  0 


or  w  =  0 


=  0  or  w  =0 


Along  <f>  -  <f>  and 


0  or  v  =  0 


(26) 
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N0x  =  0  °r 

u  -  0 

(27) 

sin  a 

" (x  +  (x  M<£,  >,«  /“  +  V  wf0/x  sin  a  +  Nx^,0  w>x  -  0 

or  w  -  0 

(28) 

m^6  -  0  or 

* 

%- 

l 

o 

(29) 

and  at  the  corners  of  the  segment 


or  w  «  0 


(30) 


4.  FORCES  AND  MOMENTS  IN  STIFFENED  CONICAL  SHELLS 

In  the  cross-  section  of  the  shell  the  strain  is  assumed  to  vary  linearly  as  in  Ref.  5, 
fx  (z*)  =  fx  -  z**x/a 

<</,  ~  *■*«<{,/* 
y*4>  <z*)  -  2z*kx<()/ a 

where  z*  is  the  physical  coordinate. 

The  stress-strain  relations  are: 

°x  (z)  =  [E/(l  -  v2)3  [<x  (z)  +  ^^(z)] 

o^(z)  =  [E/(l  -  v2)]  [(<f>  (z)  +  vrx(z)] 

rx<t,(z) =  (£/2(l  +  y*^(z) 


(32) 
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if  the  stresses  perpendicular  to  the  surface  of  the  shell  are  neglected. 

For  the  stiffeners  the  following  assumptions  are  made: 

1.  The  stiffeners  are  ,, distributed”  over  the  whole  surface  of  the  shell. 

2.  The  normal  strains  (z)  and  t ^(z)  vary  linearly  also  in  the  stiffener  according  to 
equations  (3  1).  The  normal  strains  in  the  stiffener  and  in  the  sheet  are  equal  at  their  point  of 
contact. 

3.  The  stiffeners  do  not  transmit  shear.  The  shear  membrane  force  is  carried  entirely  by 

the  sheet. 

4.  The  torsional  rigidity  of  the  stiffener  cross-section  is  added  to  that  of  the  sheet.  (The  actual 
increase  in  torsional  rigidity  is  larger  than  that  assumed.) 

Substitution  of  Eqs,  (31)  into  (32),  and  taking  into  account  the  above  assumption,  yields  the  follow¬ 
ing  expressions  for  the  normal  stresses; 
in  the  sheet 

°  x  (z)  =  [E/(l  -  v2)]  [(%  -  z**x/a  +  v  (r^  -  z*  *0  /  a)] 

a4>  (z)  =  tE /(!  ~  v2)]  [ ( ^  -  z*«^/a  +  v  i(x  -  z*Kx/a)]  (33) 

and  in  the  stiffeners 

ox(z)  =E1<,(z)  =  E1(<i-z*/(j/a) 

<V(Z)  =  E2  V(Z)  =  E2  (V  zW/a)  (34) 

where  Ej  and  F-2  are  the  "effective”  moduli  of  elasticity  of  stringers  and  frames,  respectively,  defined 
in  Section  17. 

The  membrane  force  per  unit  length  is: 

+h/2 

Ni  =  /ffx  d  z*  =  f  oxdz*  +  (l/b0x)  f  axd 

— h/2  Aj 

+h/2 

=  /  [E/(l  -  v2)][(fx-r  z* Kx/a)  +  v  ((j  -  z*k j/a)]  +  (l/bQx)  /E^  -  z*xx/a)dA1 

-h/2  A 
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*  [Eh/(  1  - „2)]  (tx  +l/e^)  +  E1  Aj  ex/bQx  -  Ejfij  (35) 

All  the  forces  and  the  moments  about  the  line  of  reference  (the  middle  surface  of  the  sheet)  are 
obtained  similarly: 

N*  =  [Eh/(1  -  v2)]  [^(l  +  zij/x)  +  ve^  -Xl  kx/x ] 

Nyb  =[Eh/(l-v2)]  [^(U^)  +  vtx  -  x2  (36) 

n*<£  =  N9b*  =  Eh  yx^/2(1  +  v) 

=  (-D/a)[(l  +  tj01  /x)kx  +  vk^  -  4  (x  /x) 

M<?S  =  (—D/a)  [(1  +  Vq2)  +  Vkx  -4 
=  (+  D/a)  [(1  -  v)  +  J?tl  /x  ] 

M<£x  =  (“  D/a)[(l  -  v)  +  ?ji2  ]kx^  (3?) 

Where  //j/x  and  /r2  are  the  increases  in  effective  cross-sectional  area  of  the  shell  due  to 
stringers  and  frames  respectively,  defined  by 

=  (1  -  v2)  Ej  Aj  /Eb0h 

ti2  =  (1  -  i/2)  E2  A2/Ea0h  (38) 

Yj/x  and  Y2  are  the  changes  in  extensional  stiffness  caused  by  the  eccentricities  of  stringers  and 
frames  respectively,  defined  by 


\1  •=  (1  —  v2)  Ej  Aj  ej  /E  bQh  a 
V2  =  (1  -  v2)  E2  A2  e2  /E  aQh  a 


(39) 
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t)Q1/x,  r)Q2  ,  r?tlA  and  r/t2  are  the  increases  in  bending  and  twisting  stiffness  of  the  shell  due  to 
stringers  and  frames  respectively,  defined  by 

*»oi  =  Ei  *01  /bo15 
*?02  =  E'2  ^02  /a0D 
'l.l  =  GlItl/b0D 

1,2  '  G2I„/a0D  H0> 

and  £  /x  and  are  die  changes  in  bending  stiffness  caused  by  the  eccentricities  of  stringers  and 
frames  respectively,  defined  by 

=  Ei  Aieia/boD 

4  =  E2A2e2a/aoD  (41) 

Nothing  has  yet  been  said  about  the  manner  of  distribution  of  the  stiffeners  and  their  magnitude. 
The  distance  between  the  frames  can  be  some  function  of  x,  and  the  magnitude  of  the  stiffeners  can  be 
a  function  of  x  and  </> .  The  changes  in  the  stiffnesses  in  Eqs.  (36)  and  (37)  are  then  functions  of  x 
and  <j>.  This  uoes  not  cause  fundamental  difficulties  if  these  functions  are  known  from  the  beginning. 

In  the  following  analysis,  the  frames  are  equal  and  equally  spaced,  and  the  stringers  are  equal. 

The  distances  between  the  stringers  are  linear  functions  of  x  (b=  bQx),  but  this  has  already  been  taken 
into  account  in  Eqs.(36)  and  (37).  Hence,  all  the  values  in  Eqs.(38)  to  (41)  are  constant. 


5.  INTERNAL  FORCES  AND  MOMENTS  AS  FUNCTIONS  OF  DISPLACEMENTS 

Substitution  of  Eqs.  (1)  and  (2)  into  Eqs.(36)  and  (37)  yields  the  additional  internal  forces  and 
moments  acting  during  buckling  as  functions  of  the  additional  displacements: 

=  [  Eh/(  1  -  v-2)]  [(1  +  //1/x)u_x  +  r-fv^/x  sin  a  +  u/x  -  wcota/x)-  yj  wxx/x] 
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=  [Eh/(1  -  i/2)]  [(1  +  fi2)(,v^/x  sin  a  +  u/x  —  w  cot  a/x)  +  vu  x  -  y2  (w  x  /x 

+  w,^/x2sin2  a)l 

Nx<£  *  [Eh/2 ( 1  +  i/)]  [v>x  -  v/x  +  u^/x  sin  a]  (42) 

M*  "  [-D/a][(l  +  ?j01/x)  w>xx  +  v(w^x/x  +  w^^/x2  ain2a)  -  CjU^/x] 

M^6  =  f-D/a]  [(1  +  7jfl2)  (w>x/x  +  w^^/x2  sin2  a)  +  v  w  xx 

“  ^2  X  S*n  a  +  u/x  —  w  cot  fl/x)] 

^x<^>  “  [+  D/a]  [(1  -  v)  +  jjtl  /x]  (w^/x  sin  a  -  w^/x2  sin  a) 

=  [~D/a]  [(1  -  v)  +  r)t2  ]  (w(X?i/x  sin  a  -  w^/x2  sin  a)  (43) 


6.  TRUNCATED  CONICAL  SHELL  UNDER  HYDROSTATIC  PRESSURE. 

A  simply  supported"  circular  truncated  conical  shell  is  considered,  closed  at  the  ends  by  bulk¬ 
heads  which  offer  no  restraint  against  displacement  or  rotation  of  the  generators  of  the  shell  while  being 
rigid  perpendicular  to  them. 

The  load  is  uniform  hydrostatic  pressure  acting  on  the  shell  and  the  bulkheads.  It  is  assumed  that 
the  stress  state  prior  to  buckling  is  represented  satisfactorily  by  the  membrane  stresses 

Nx0  =  -(p/2) ax  tan  a 

N9io  =  -P  ax  tan  a 


N^0=<> 


(44) 
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In  the  complete  truncated  conical  shell  the  boundary  conditions  along  the  lines  cf>1  =0  and 
^2  “  2n  (which  are  the  same  line)  are  satisfied  automatically,  and  the  shell  has  no  corners.  For  the 
case  of  hydrostatic  pressure,  Eq.  (18)  becomes  then 

2ff  r2 

SV  =  -  /  f  *[(xNx).x/ax  -  N0/ax  +  Nx^,x/ai  sin  a]  Su+  [N^/ax  sin  a 

+  +  [(xMx)fiI/a2x  -  M^x/a2x  +  M^/ax2  sin2  a 

~  (XV.x</a2*2  sin  a  +  (X  *2  3*n  °  +  COt  “/aX 

-  p  tan  a  (x  wxi/2  +  w^x  +  w^^/x  sin2  a)]  8w)  a3  x  sin  a  dxd<£ 

2ff 

+  /  !ax[Nx  aSu  +  Nx^aSv-Mx5(Wfi)]  +  [(xMx)fi  _  -  M^/sin  a 

X5=*2 

+  M0x,0/sin  a  -(p/2)a2x2  tan  a  wj  a  8  w  1  sin  a  d<f>  =  0  (45) 

x=l 

The  critical  pressure  is  obtained  from  the  above  condition  that  the  virtual  work  must  vanish. 

After  substitution  of  Eqs.  (42)  and  (43),  Eq.  (45)  becomes  a  function  of  the  displacements  only.  Thus, 
the  problem  reduces  to  finding  displacements  u,  v,  and  w  which  satisfy  Eq.  (45). 

7.  DISPLACEMENTS. 

The  admissible  displacements  used  in  the  solution  of  Eq.  (45)  are  assumed  as  in  Refs.  1,  6  and 
7  in  the  form 

N 

u  =  !Jim  2  A  x®  sin  t<f> 

n  =  l 
N 

v  «  $Sm  £  Bn  x®  cos  t^ 


IS 


N 

w  =  5m  £  C  x“  sin  t<£  (46) 

n  =  l 

Where  5m  indicates  the  imaginary  part  of  the  expressions,  An  and  Bn  are  complex  coefficients 
and  C  and  t  are  real,  (t  is  the  number  of  circumferential  waves  of  the  buckling  deformation),  s  is 

n  v 

the  complex  number 

s  t»  y  +  in  (47) 

y  and  are  real,  and  are  calculated  from  the  boundary  conditions,  n  is  a  real  integer,  and 

i=V(~D. 

The  displacements  of  Eqs.  (46)  solve  the  first  two  stability  equations  of  an  unstiffened  conical 
shell  exactly,  and  the  third  stability  equation  is  solved  in  Ref.  1  by  the  Galerkin  method. 

In  the  case  of  a  stiffened  conical  shell  these  displacements  do  not  solve  the  first  two  stability 
equations  exactly.  Nevertheless,  since  it  is  assumed  that  the  form  of  buckling  in  general  instability  of 
a  stiffened  conical  shell  does  not  differ  much  from  that  of  an  unstiffened  conical  shell,  the  displacements 
of  Eqs.  (46)  are  used  also  here. 


8.  EVALUATION  OF  /9  AND  y  BY  COMPLIANCE  WITH  BOUNDARY  CONDITIONS 


The  shell  is  assumed  to  be  simply  supported.  At  the  end  planes  (x  =  1,  x2),  the  radial  displace 
ment  w  and  the  longitudinal  moment  Mx  must  vanish. 

A  typical  term  of  the  displacement  series  (46)  is 


w 


n 


y+in/9 

5m  Cn  x®  sin  t^  =  5m  Cn  x  sin  t 

y  t 

5m  Cn  x  [cos  (n  /3  In  x)  +  i  sin  (n  /3  In  x)]  sin  t<£ 


(48) 


Hence,  its  imaginary  part  is 

y  , 

wb  =  Cn  x  sin  (n  fi  In  x)  sin  t  cjy 


(49) 
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Since  In  1  =  0,  w  vanishes  identically  at  the  small  end  plane  (x  =  1)  of  the  cone.  In  order 
that  w  should  also  vanish  at  the  large  end  plane  (x  =  x2),  p  must  be  defined  as 


P  =  jr/ln  x2  (50) 

From  Eqs.  (43)  the  longitudinal  moment 

=  -  (D/a)  [(1  +J?01A)  w  xx  +  v(w  x/x  +  w  ^/x2  sin2  a)  -  u>x/x  ]  (51) 

At  the  boundaries  considered,  w  is  zero  and  u  x  is  very  small  and  may  therefore  be 
neglected.  Hence,  at  x  =  1,  x„  , 

(1  +  rjQl/x)  w  xx  +  v  w>x/x  =  0  (52) 

Substitution  of  Eq.  (49)  into  Eq.  (52)  yields  for  the  typical  displacement  term 

x  2  sin  (n/91n  x)  |  [y  (y  -  1)  -  n2  $  ]  (1  +  »?01A)  +  vy  I 

(53) 

y— 2 

+  x  cos  (n  p  In  x)  n  /9  [(2y  -  1)  (1  +  )?01A)  +  v]  ==  0 


and  after  substitution  of  the  boundary  conditions 

(2y-  1)(1  +t)01/x)  +  v  =  0  (54) 

Simultaneous  fulfilment  of  condition  (54)  at  both  boundaries  is  impossible.  It  can  be  fulfilled 
only  in  one  of  the  boundaries,  say  at  x  =  1,  where  Eq.  (54)  becomes 

(2 y  -  1)  (1  +  r?01)  +  v  =  0  (55) 

Hence 


y  =  (1/2)  [1  -  v/(l  +  »701)] 


(56) 
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Equation  (56)  fulfils  the  boundary  condition  =  0  only  approximately  for  stiffened  shells.  In 
shells  stiffened  by  stringers,  small  internal  longitudinal  moments  act  at  the  boundaries.  However,  in 
shells  stiffened  only  by  frames,  the  equilibrium  boundary  conditions  of  zero  longitudinal  moments  is 
satisfied  rigorously.  Since  in  the  absence  of  stringers 

i=0 

(57) 

=0 

and  Eq.  (51)  yields,  therefore, 

w  +  v  w  /x  =  0 

,11  ,z  ' 

By  substitution  of  Eq.  (49)  into  Eq.  (58)  one  obtains  then  for  both  x  =  1  and  x  =  x2 
(2y  -  1)  +  */  =  0 

and  hence 

y  =  (1  -  v)/2  (60) 

The  same  result  could  be  obtained  by  substitution  of  rjQ  x  =  0  into  equation  (56).  The  value  of  y  for 
conical  shells  stiffened  by  frames  only,  is  the  same  as  that  for  unstiffened  conical  shells  (Ref.  1),  be¬ 
cause  the  frames  do  not  appear  in  the  expression  of  the  longitudinal  moment.  Hence  Eq.  (60)  fulfils 
the  equilibrium  boundary  condition  of  zero  longitudinal  moments  for  unstiffened  and  frame-stiffened 
conical  shells. 


(58) 


(59) 
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9.  MODIFICATION  OF  THE  FIRST  TWO  STABILITY  EQUATIONS 

At  the  boundaries  w  and  S  w  vanish.  Hence  also  the  boundary  term  multiplied  by  8  w  in 
Eq.  (45)  is  zero.  Substitution  of  Eqs.  (42)  and  (43)  into  the  first  two  stability  equations  of  Eq.  (45) 
yields : 

2/r  *2  , 

SU  =  -  f  f  |[Eh/a  (l-i/2)]  |u  xx  (l  +  /i,/x)  +  u  x/x-u(l  +  *t2)/x2+(l-v)  u>(^/2x2  sin2  a 
0  1 

+  (1  +  v)  v  x<^/ 2x  sin  a  -  [(3  -y)/2  +  n2]  v^/x2  s‘n  a~Xi  W/XXXA 
+  x2  w  x/x2  +  x2  w^^/x3  sin2  a  -  v  cot  a  w  x/x  +  (l  +  p2)  cot  a  w/x2|  Su 

+  [Eh/a  (1-  v2)]  I  (1  +  u)  u  x^/ 2x  sin  a  +  t (3  —  v)/2  +  u^/x2  sin  a 

+[(1  -  v)/2]v>xx  +  (1  -  v)  v_x/2 x  +  (1  +  n2)  V,<^A2  sin2  a 

-  (1  -  v)  v/2x2  -  y2  w^^/x3  sin3  a  -  x2  w^A2  sin  « 

-  (1  +  n„)  cot  a  w  ,/x2  sin  a  1  8 v  +  [(xM  )  ../a2  x 

-  M^.x/a2  X  +  V^/a2  x2  sin2  8  -  *2  8in 

+  (xM^x)  x(^/a2  x2  sin  a  +  cot  a  N^/ax  -  p  tan  a  (x  w  xx/2  +  w  x 

t  w  s*n2  &3  x  s*n  0  ^x  ^ 

2rr  *=x2 

+  f  {ax  [Nx  aSu  +  Nx<^  a  S  v  -  8  ( w  x)]j  sin  a  d<fi  =  0  (61) 

The  displacements  of  Eqs.  (46)  do  not  solve  the  first  two  stability  equations  exactly.  In  order 
to  “correct”  this,  a  set  of  terms,  whose  sum  is  zero,  is  added  to  Eq.  (61).  Then 

2!T  x  , 

SU  =  ~  /  f  3  [Eh/a  (l-i,2)]|(l+Ml/x)u  xx+u  x/x-(l+^2)u/x2+(l-v)u((^/2x2  sin2  a 
0  1 
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+  (1  +  v)  x  sln  a~  t(3-v)/2+//2]  v^/x2  sin  a-xx  w>xxxA  +  X2  W(XA2 

+  x2  w  (^/x3  sin2  a  -  v  cot  a  w  x/x  +  (1  +  n2)  cot  a  w/x2j  Su 

2  2 

+  [ E h/a  (l-^2)][Ml  n>xxA1-/t1  u.xx/ki-VixWfXXXA2+X1xWfX„/k2 

+  X2  w.x/xk3  -X2  W.x/X  k3  +X2  W,^/k4  sin2  «  x"  "X2  W,^/k4  sin2  ax2]  Su 
+  [  Eh/a(l  -i/2)  ]  1  (1  +  v)  u  x(^/2x  sina  +  [  (3~r/)/2  +  n2]u^/x2  sin  a  +  (l-v)  v(XX/2 

+  (1  -  v)  v  x/2x  +  (1  +  /i2)  v^/x2  sin2  a  -  ( 1  -  v)  v/2x2  - x2  W.^/x3  sin3  a 

-  x2  w  ^/x2  sin  a  -  (1  +  /r2)  cot  a  w  ^/x2  sin  a  !  Sv 

+  f  E  h/a  (l-i/2)  ]  [-  y2  w^^</(/k5  x 2 sin3 a  +  y2  w>(, AB  * 2  sin3  a 

-  *2  w.x/k«  x  sin  «  +  X2  w,x</k6  x  sin  a]  5v 

+  t(xMx),„/a2x-M^.x/a2  X  +  Mvi,,^/a2  x2  sin2  W’  x2  sin  a 

+  (x  M^)  ^/a2  x2  sin  a  +  cot  aN^/at -p  tan  a  (x  w  xx/2  +  w  x 
+  w  ^^/x  sin2  a)]  8w  j  a3  x  sin  a  dx  d</> 


277 


+  /  |  ax  [N  aSu  +  N,aSv  —  M  5  (w  )]  |  sin  a  d<f>  =  0 

J  x  X  G>  X  '  ,x'  . 

0 


(62) 


The  expressions  multiplied  by  the  second  Su  and  Sv  are  these  new  terms.  After  rearrangement  of 
the  terms  of  Eq.  (62)  one  obtains  : 


2  77  x 


SU 


o  t 


■{ 


=  -/  /  |[Eh/a(l-l72)]l(l+,r1A1)UiXX  +  UiX/x-(l+/i2)u/x2+(l-v)u^/2x2sin2a 


+  (1  +  v)  V  x<^,/2X  sin  a  - [(3  -  u)/2  +  fi2]  v^/x2  sina-Xj  xwfXXXA2  +X2  W,XA3  x 

+  x2  w,00/  ^4  x2  s'n2  a-vcotaW(X/x  +  (l  +  fr2)cotaw/x2  1  Su 
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+  [  Eh/a ( 1  -  „2)  ]  [„,  «„/*-/»!  w,,„/*  +  Vi  *  »„„aJ 

+  V2  w,x/*2  "  V2  w,, Asx  +  *2  w (^/x*  sin2  «~X2  Wf<^/k4  x2  sin*  a]Su 
+  [Eh/a(l  —  w2)]  ((l  +  v)uiX<^/2*  8ina  +  [(3-v)/2  +  ^2]u  ^/x2sin«+(l-»/)v  %  J2 
+  ( 1  - »/)  v  x/2x  +  (1  +/ij)  V,^A2  sin2  a- (1  -  v)  v/2x2  -  V2  w^^/kjX2  sin*  a 
~  X2  x  sin  a  -  (1  +  /r2)cot  a  w  ^/x2  sin  a]Sv 

+  [Eh/a(l-v2)][-V2W>^/x*  8in3a  +  x2w(^0/k5x2  sin*  a 
“X«  W,x^/x2  8in  a  +  v2  w,x0/k«  X  8in  a]  Sv 
+  [(x  M*\„/a2x  -  M^,x/a2x  +  M^,^/a2  *2  ain2a  -  (xMx0\*ya*  x2  8in  a 

+  ^xM^*\x</(/ ®2x2  sin  a  +  cot  a  N^/ax  -  ptan  a(x w>xx/2  +  w 
+  w,^  /  x  sin2  cn)  ]  5  w|  a3  x  sin  a  d  x  d 

2ff  x~x2 

+  /UxIN^iu  +  NX<^  afiv  -  Mr5(w  x)]  |  sinad^-O  (6 

0  x=l 


Now,  in  Eq.  (63)  the  expression  multiplied  by  the  first  Su  is  the  first  “corrected  stability 
equation"  and  the  expression  multiplied  by  the  second  Su  is  the  “error-load”  of  the  first  stability 
equation.  In  the  same  way,  the  expression  multiplied  by  the  first  Sv  is  the  second  “corrected  stabi¬ 
lity  equation"  and  the  expression  multiplied  by  the  second  Sv  is  its  “error-load",  kj.kj.kj  ,k4,kg  and 
kfi  are  the  correcting  coefficients",  which  will  be  calculated  by  equating  the  virtual  work  done  by  the 
“error-loads”  to  zero. 

Ks  a  matter  of  fact,  for  calculation  of  the  “correcting  coefficients*  it  is  usually  enough  to  con¬ 
sider  only  the  first  term  of  the  displacements  series. 
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10.  SOLUTION  OF  THE  FIRST  TWO  “CORRECTED”  STABILITY  EQUATIONS 

The  following  notation  is  introduced  for  brevity: 

Ft  /kj  =  /iji 

2 

^•1^2  “  Yj2  Y2//^S  =  ^25 

Y2A3  =  X23  ^2^6  =  *26 

X2/k4  =  *2  4  (64) 

Since  the  displacements  of  Eqs.  (46)  solve  the  first  two  “corrected”  stability  equations  of  Eq.(63) 
exactly,  they  may  now  be  written  for  these  displacements  as 

(l  +  Fn)u|XX  +u,*/*  -  (l  +  <x2)u/x2  +  (1  -^u^/^x2  sin2  a  +  (1  +  v)v  x^/2x sin  a 

-  [(3  -  »')/2  +  /i2] v^/x2  sin  a-x12^wixxx  +  k23wix/x  +  x24wi(/)yx2ain2a 

-  v  cot  a  w  /x  +  (1  +  /ij)  cot  0  w/12  =  0 

(1  +  v)  uiX(^/2  x  sin  a  +  [(3  -  y)  /2  +  q2]u  ^  /x2  sin  a  +  (l  -v)v  fXX/2  +  (\-  v)vx/2x 
+  (1+  sin2  a  -  (1-v)  v/2  x2  -  v/^^/x2  sin3  a 

“  *26  W,x^/;|C  sin  a  -  (1  +  H2)  cot  a  w,4>/'x2  sin  0  ”  0  (65) 

and  after  substitution  of  the  displacements  of  Eqs.  (46)  into  them  one  obtains 

w  r 

2  xs  2  sin  t  3  A  [(l  +  #/  )  s  (s  -  1)  +  s  -  (1  -  i/)t2  /2  sin2  q] 
n— 1  L  11 
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+  Bn  j  -  (1+  v)  t  a  /2  sin  a  +  [(3  -  v)  /  2  +  n2  ]  t/  sin  a  } 

+  Cn  \i2  a(s-l)(a-2)+Xa,  8  “  *24  t2/sin2  a 
-  v  s  cot  a  +  ( 1  +  n2)  cot  a  ]  ^  =  0 
N  , 

2  x*- 2  cos  t  cf>  3  An  i  (1  +  v)  st/2  sin  a  +  [(3  -  v) /  2  +  ji2]t/sin  a  i 

n  =  l  *■ 

+  Bn  [(1-v)  s(s  —  l)/2  +  (1-v)  s/2  -  (l  +  ^2)t2/sin2  a  -  (l-v)/2] 

+  Cn  [  t3/sin3  a  -  y26  t s/sin  a  -  (l  +  ^2)t  cot  a/sin  a]|  =  0  (66) 

For  Eqs.  (66)  to  be  satisfied,  every  term  of  their  series  must  be  zero.  After  division  by  the 
factors  x“-2sint<£  and  x*-2  cost<ji,  which  multiply  every  term  of  the  series  in  the  first  or  the 
second  equation,  respectively,  algebraic  equations  are  obtained.  The  solution  of  these  algebraic 
equations  yields  the  coefficients  An  and  Bq  as  functions  of  Cn 

A  =  (A  +  i  A.  )C 

n  rn  in'  n 

B  -  (B  +  iB.  )C 

11  '  rn  in'  n 

where  A  .  A,  ,  B  ,  B.  and  C  are  real  coefficients  defined  by:  - 

rn 1  in'  rn'  in  n 

A„  =  (MqM1  +  N0N1)/(m"  +  N02) 

Ai„  =  <MoNi  ~NoMi >/(Mo  + 

Br„  -  <M0  M2  +  ^0  N2>/(M0  +  Nc> 

Bin  =  (M0N2  -N0M2)/(M2  +  N20) 


(68) 
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where 

M0  (n)  =  M01  t4  +  M02  t2  +  M03 
N0(n)  -  Nfll  t2  +  Nq2 

M,  (n)  =  Mu  t4  +  M 12  t2  +  Mj3 

N,  (n)  =  N„  t4  +  N12  t2  +  Nj3 
M2  (n)  =  M21  t5  +  M22  t3  +  t 

N2  (n)  -  N21  t3  +  N22  t  (69) 

Moi  (n)  =  doi  /ain4  a 

M02  (n)  =  ('/sin2  a)  [n2  ft*  (~dQ2)  +  y2  dQ2  +  y  +  dM  ] 

M03  (n)  =  n4  ^  d05  +  "2  P2  (-6  y2  d05  -  3  y  d06  -  d07>  +  ^  d05  +  V®  d06 

+  V2d07  +  Vd08  +  d09)  (70) 

N01  (n)  =  (n/3/sin2  a)  (2y  dQ2  +  dQ3) 

Nq2  (n)  =  n3/)3(-4yd05  -  dog)  +  n/3  (4y3  dQ5  +  3  y2  dQ(5  +  2y  dQ7  +  dQ8)  (71) 

Mn  (n)  =  (1/sin4  a)  (y  du  +  dJ2) 

M12  (n)  =  (1/sin2  a)  [n2  fi2  (-3  y  d13  -  du)  +  (y3  dJ3  +  y2  du  t-  y  d,5  +  dlfi)l 

M1S  (n)  =  n4  /34  (5  y  dj  ?+  dt  g)  +  n2  ft*  (.-  10  y3  d17-  6y2d,B  -  3yd19  -  dno) 


+  (y5  d1?  +  y4  d18  +  y3  d19  +  y2  dUfl  +  ydin  +  d,12) 


(72) 
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Nu  (n)  =  n  f3  dn  /sin4  a 

N12  (n)  =  (1/sin2  a)  [n3 /33(-d13)  +  n£(3  y2  d13  +  2  y  dM  +  dls)] 

N13  («)  =  «5/S5d17  +  n3/33(-10y2d17  -4yd,,  -d,,) 

+  n  /3(5y4d17  +  4y3d18+  3y2di9+2ydiio  +  din)  (73) 

MJt  (n)  =  d21  /  sin*  a 

^22  (n)  “  O/sin3  a)  [n2  ^  (-d22  )  +  (y2  d^  +  y  d23  +  d^)] 

M23  (n)  “  (1/sin  a)  [n4  /3  d^  +  n2  /S  (-6  y2  d&  -  3  y  d2ft  -  d2? ) 

+  V  d25  +  y3  +  y2  d27  +  y  d28  +  d29>l  (74) 

N21  (n)  =  (n  /3/sin3  a)  (2  y  d^  +  d^  ) 

N2i2  (n)  =  (1/sin  a)  [n3  03  (-4  y  d^  -  d26 )  +  n  0  (4  y3  d^  +  3  y2  d^ 

+  2  y  d27  +  d2«)3  (75) 

doi  =  0  “  v)  (1  +  //2) 

d02  =  2  v  —  2(1  +  n2)  (1  +  /tu) 

d03  ~  ^  ^11  ^1  +  Mj) 
dw  -  -2(1  -  v)  (1  +  n2) 
d05  -  (l-v)(l  +  „11) 
d06  =  ~  (1  ~  ^  Mu 
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d07  =  -(1  -  ^  (^11  +  //2  +  2) 
d08  =  (1  -  ^11 
d09  =  (1  -  ^  d  +  H2) 

dn  =  -  d  +  *")  *25 

di2  =  -  (1  +  //2)2  +  (3  -  v  +  2 n2)\x 

d13  =  -2(1  +  /^2)  Via 

dU  =  6d  +  ^2^12  +  d  -  v)X2t  +  ^  +  ^^26 

d15  +  /i2)  X'^  +  2^l  +  M2)X23-^_>/  +  2^2^  ^26  +  (1  ~  v)  (1  +  /i2)cot  a 

di6  =  ~  v)  x2i  -  (1  -  v)  (1  +  ft2)  cot  a 

d  17  =  (1  “  ^Vl2 

d18  =  -3(l-„)X'la 

di9  =  ^  ~  ^  X'  12  -  (1  -  ^)  x'23  +  t'(l-i')  cot  a 
dUo  =  ^  '  ^Xi2  “  d-  v)  U  +  n2)  cot  a 
diu  -2(1  - v) \ 12  4  (l-t/lx^  ~  I'd  -  v)  cot  a 
dil2  1=1  (1  -  p)  (1  +  /j.2)cot  a 

d21  =  (1“^X25 
d22=-2(l+,u)V25 
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d23 

=  2^11 

*25 

(1  + 

£ 

+- 

1 

M 

1 

* 

*26 

d24 

=  -(3- 

*'  +  2  H2 

^*24 

+  2(1+  n2)v25 

t  (1-v)  (I  +  n2)  cot  a 

d25 

=  -(1  + 

v)  vla 

d26 

=  2(2u 

-  M2>  *12  + 

2  *26  U  +  f*n> 

d27 

=  -2/in 

*26  + 

(7- 

5»/+  6^2)Vl2  + 

(1  +  ^  V23 

+  [ 2(1  +  n2)  (1  +  ^u)  _  \/(\  +  »/)]  cot  a 

d28=~  ^  +  ^2^  cot  a  ~  2  (3  —  v  +  2/i2)  ^12  +  (3  -  v  +  2 li2)\23 

-  2(1  +  ^2)  x26  +  (1  -  v)  (1  -  v  +  n2)  cot  a 

d29  =  (1  -  u)  (1  +  n2)  cot  a  (7g) 

Now,  since 

x*  =  x>'+1'’  0  =  xV[cos  (n/3  lnx)  +  1  sin  (n /9  In  x)]  (79) 

the  displacements  become,  after  substitution  of  Eqs.  (67)  into  Eqs.  (46) 

N 

u  =  3m  £  Co  (Arn  +i  Aia)  x^  (cos  (n  /3  In  x)  +  i  sin  (n  f}  In  x)]  sin  t 

n=l 

N 

V  =  Sm  £  Cn(Btn  +  i  Bin)xy  [cos  (n/3  lnx)  +  i  sin (n/9 lnx)] cost <f> 

n — I 

w  =  $m  £  ( cos  (n  /3  In  x)  +  i  sin  (n  /3  In  x)  ]  sin  t  d>  (80) 


The  imaginary  part  of  these  equations  is 
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N 

u  =  sin  t $  2  xy  [Afn  sin  (n  /S  in  x)  +  AJn  cos  (n  /9  In  x)]  Cn 

n  — 1 

N 

v  =  cos  t <f>  2  xY  [Brn  sin  (n/91n  x)  +  B(n  cos  (n/31n  x)]Cn 

n=l 

N 

w  =  sin  t<£  2  x.y  [  sin  (n  fi  In  x)]  Cn  (81) 

n  — 1 

where  all  the  values  are  real. 

The  three  displacements  are  thus  expressed  as  functions  of  one  set  of  real  arbitrary  coefficients 
Cn.  Hence  the  variations  of  the  displacements  are 

N 

2u  =  sin  2  xV  [Afn  sin  (n/9lnx)  +  A(n  cos  (n/91n  x)]  5Cn 

n=l 

N 

Sv  =  cos  t(p  2  xY  [Brn  sin  (njQlnx)  +  BJn  cos  (n/31nx)]  SCa 

n  =  l 
N 

Sw  =  sin  t<£  2  xY  [sin  (n/3  In  x)]  SCB  (82) 

n  =  l 


11.  CORRECTING  COEFFICIENTS 


In  Section  9  “correcting  coefficients”  were  introduced  which  then,  in  Section  10,  permitted  the 
expression  of  the  u  and  v  displacements  as  functions  of  the  radial  one,  w.  These  “correcting  coef¬ 
ficients"  will  now  be  evaluated  by  equating  the  virtual  work  done  by  the  error-loads  of  the  first  two  sta¬ 
bility  equations  of  F.q,  (63),  to  zero. 


2jt  x 


0=-  /  /l  [Fh/a  (1—  v2)]  [u(ix  H/x  -  UfiiMl/kj  -  Xx  /x  +  Xl  x  W(MJ| /v\ 


n  1 


+  V2  W,*/x2  ~  ^2W,>/k3  x  +  X'2W,«A<A/xS  Sin2a~^2W,0d./k4  x2  8in  a]  5U 
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+  [Eh/a  (1  —  v2)]  [-XjW^^/x3  sin3  a  +  Y2 x2  sin3  a 
-  X2  wlXj,/x2  sin  a  +  y2  WiX^/kfi  x  sin  a]  S  v  |  a3  x  sin  a  d  x  d  <p  (83) 


If  all  corresponding  pairs  are  equal  to  zero,  the  entire  Eq.  (83)  vanishes  too.  The  "correcting 
coefficients”  may  be  calculated  in  this  manner.  For  example,  for  one  obtains 

2(T«2 

-  /  /  [Eh/a  (1  -  i/2)]  [ujXx  til  /x  -  U)ix  til  /kj  ]u  a3  x  sin  a  dx  d0  =  0  (84) 

o  1 

and  then 


f  f  u  xx  u  x  dxd<£ 
o  1 

k1= -  (85) 

277  x2 

/  /  uxx  udxd0 
0  1 


Similarly  the  remaining  “correcting  coefficients”  are  obtained 


2nxl 

f  f  w  xxx  u  x2  dxd  i/> 


f  f  wrxixudxd9i 
0  1 

2JTX 

/  j  w  u  dxd(£ 
o  1 

277  *2 

f  f  (WfX/x)udzd0 
0  1 
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/  /  (w.^/*2)  dxd  <£ 


/  f  (w,^/x)vdxd«^ 


f  f  (w.4dd>/x2)v  dxd</> 


/  /  w^vdxd^ 


/  f  (w,/*)  v  dxd  ^ 
0  1 


Substitution  of  the  displacements,  Eqs.  (81),  into  Eqs.  (85)  and  (86)  yields: 


2,  S_.C“C“tIiy  1(n-m)AcoA-+ir  1  (rain) Ac „  Aim  + I***  1(»,m)AdBAra+i;ri(n,m)AiBA1J 


n=  1  m=l 


s  1  2MAclAB+i;y  2(m,n)AcnAim  +l22y  2(n,m)AdnArm+I32y  %,m)A.  A.l 


S  S  Cn  Cm  (n>m)a3n  Arm  +  (m  -n)  a3  n  Aim  +  ^  (“>m)  b3  n  Atm  +  ^  '(n<m)b A 


2  n=l  m  =  l 


n^^=iCaC,"[Il  (m’n)a3nAi,»  +  ^  “  (»•*) b3n  Atnl  +  l\*  ' 
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N  N 


S  S  CnCnJ1iy  1(n.m)yArD1  +  I2r'1(m'n)>'Alm  +  I!y_1(n'm)n^Arm  +  Isy  1  («».">) n  0  A,J 


n  =  l  m  =  l 


(89) 


N  N 


£  2  Cn  Cm  [I^^~i(n,m)y  Arm  +  l*r  \m,a)y  Aim  +  l‘2'  ‘  (n,m)n  0  A^'  ‘(D.mJn^AJ 
n*=l  m  =  l 


r2r-2, 


t2y-2 


T2y-2 


N  N 


S  S  cn  +  f‘MAim] 


k4  = 


n  =  l m=l 


S  S  Cn  Cm  [if  2  (n,m)  Arm  +  1^*  (m-n)  Aim  1 

n=l m— 1 


(90) 


k,  = 


2  2  CnCm  [if -1(n,m)Brm  +  if  WjB.J  1 

n  =  l  m—  1 


N  N 


S  2  CnCm[lf"2(n,m)Brni+I2y"2(m,n)Biin}  af"* 


n=l m=l 


(91) 


N  N 
2  2 
n=l m— 1 


CnCm[lf-1(n,m)yB| 


ra.+I22y~1(m>n>yBi 


im+  ir’(“.nl)nP[m+  Ij^V.mXl^B^ 


2  2  CnCm[lf  2(n,m)yBrni+  I^m.nVB^+lf"2  (n.m^^lf^^.m^B^ 


2y-l 


K 


6 


(92) 


2y-2 


It  is  possible  to  combine  several  correcting  coefficients  in  order  to  reduce  the  number  of  coef¬ 
ficients.  For  example,  the  coefficient  kff  was  found,  during  the  calculations,  to  be  very  small. 
Hence,  it  was  combined  with  k5  and  redefined  as 


k5  =  k6 


5  6 


(93) 


From  Eq.  (83),  one  obtains  then 
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2n  * 

f  J  a‘nS  fl  +  x  8*n  a)  v  *  dxd  <£ 


^56  " 


l  {  (W,^/X$  8inS  fl  + 


WiX0/x28ina)  vxdzd^ 


and  finally 


(t2  / sin2  a)  k2y  1  -  k2y  1 


(ta/eia*  a)  -  k‘>“! 

o  6 


The  various  terms  which  appear  in  Eqs.  (86)  to  (94)  are  defined  as  follows 


Ac  “  My-  W  -  n2  ^2]  Afn  -  n  /8  (2y-  1)  Aln 
Adn  -  My-  D  -  n2  /92)  Aln  +  n  /9  (2y  -  1)A 


a2n  =  y(y-  1)  -  n2  /9 


b2»  -  n/9(2y-l) 


a3»  =y(y-l)(y-2)-n2/92(3y-3) 
b3n  ■=  »i9(3y2  r  6y  +  2  -  n2  /32) 

K  *2 

Ij  (n,  m)  -  fxK  sin  (n  /9  In  x)  sin  (m  $  In  x)  d  x 


(1/2)  (k+1)  [1  -x"+1(-l)“+m][ 


(n  +  m)2  /32  +  U  +  l)2  ~  (n-m)2/92+U+l)2 
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K  2 

12  (n,m)  =  f  xK  cos  (n  /3  In  x)  sin  (m  0  In  x)  d  x 

1 

,  k+i  ,  ,n+m,  ,  n  +  m  m— n 

=  (1/2)  0[l  -  x  (-1)  ][ - +  - ] 

(n  +  m)2  02+(x  +  l)2  (m-n)2 02  +  U+  l)2 

K  x2 

13  (n,m)  =  f  xK  cos  (n  0  lnx)  cos  (m  0  lnx)  dx 

l 

=  (l/2)(K+l)[x*+1(-l)D+,n-l][ - - -  +  - - - 1  (98) 

(n  +  m)2  02  +  (k+  l)2  (n-m)202  +(k+1)2 

The  physical  interpretation  of  the  “correcting  coefficients*  is  that  they  introduce  artificial  ties 
between  the  displacements  u,  v  and  w,  which  reduce  the  degrees  of  freedom  of  the  shell.  The  shell 
is  thus  artificially  stiffened  and  the  critical  load  obtained  is  higher  than  the  actual  one.  For  calcula¬ 
tion  of  the  correcting  coefficients,  the  first  term  of  the  displacement  series,  which  represents  the  domi¬ 
nant  component  of  the  buckling  mode  under  hydrostatic  pressure,  is  usually  sufficient.  If  the  first  ap¬ 
proximation  yields  the  required  accuracy,  the  correcting  coefficients  have  no  effect  on  the  critical  load. 
Since  addition  of  the  second,  third  and  higher  terms  of  the  displacement  series  only  improves  the  one- 
term  solution,  the  increase  of  the  critical  load  due  to  the  correcting  coefficients  must  be  negligible. 

It  should  be  pointed  out  that  the  coefficients  A  ,  A  B  ,  and  B.  which  appear  in  the  right 
side  of  Eqs.  (87)  to  (92)  are  themselves  functions  of  the  correcting  coefficients.  Hence  an  iteration 
procedure  is  required.  One  first  assumes  values  for  the  correcting  coefficients,  and  after  calculation 
of  the  critical  load  checks  if  they  were  assumed  properly.  If  large  differences  are  obtained,  the  critical 
load  has  to  be  recalculated  with  the  new  correcting  coefficients.  Since  the  requirement  of  zero  virtual 
work  done  by  the  error-loads  of  the  two  stability  equations  is  not  a  mandatory  requirement,  but  is  only 
a  means  for  improvement  of  the  accuracy  of  the  solution,  the  correcting  coefficients  need  not  be  deter¬ 
mined  very  accurately.  The  correcting  coefficients  can,  therefore,  be  calculated  with  the  coefficients 
Arn,  A.  ,  Brn  and  B.n  obtained  from  the  approximate  solution  given  in  Section  IS,  where  the  latter 
coefficients  are  not  functions  of  the  former. 

Although  the  use  of  the  correcting  coefficients  introduces  artificial  ties  between  the  displacements 
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u,  v  and  w,  these  ties  are  not  far  from  the  actual  ones.  In  Ref.  13,  the  authors  have  shown  that  the 
rigorous  solution  of  the  Donnell  type  stability  equations  for  an  unstiffened  cylindrical  shell  under 
hydrostatic  pressure  is  valid  also  for  those  of  a  stiffened  cylindrical  shell  failing  by  general  instability, 
since  the  behaviour  of  a  stiffened  cylindrical  shell  is  similar  to  that  of  an  unstiffened  one.  It  is  reason¬ 
able  to  assume  that  the  same  will  occur  for  conical  shells.  Hence  one  may  expect  that  the  buckling  dis¬ 
placements  of  unstiffened  conical  shells  will  be  suitable  for  the  case  of  general  instability  of  stiffened 
conical  shells.  However,  though  these  displacements  solved  the  first  two  stability  equations 
exactly  in  the  case  of  unstiffened  conical  shells,  they  do  so  in  the  case  of  stiffened  conical  shells  only 
after  the  first  two  stability  equations  have  been  "corrected”  slightly  by  the  correcting  coefficients. 


12.  BOUNDARY  CONDITIONS 


The  buckling  displacements  have  to  satisfy  the  following  boundary  conditions  : 

a.  Geometrical  boundary  conditions 

w  *  0 

at  x  «=  1  and  x  =  x2 

v  =  0 

b.  Equilibrium  boundary  conditions 

Ns  =  0 

at  x  =  1  and  x  =  x, 

M  =  0 

X 


(99) 


(100) 


The  displacements  of  Eqs,  (46)  do  not  fulfil  all  the  boundary  conditions.  In  Section  8  the  values 
of  y  and  £  were  determined  from  the  compliance  with  the  boundary  condition  of  zero  radial  displace¬ 
ment,  and  also  of  zero  longitudinal  moment,  when  the  shell  is  stiffened  by  rings  only.  The  displacements 
do  not  fulfil  the  boundary  conditions  of  zero  circumferential  displacement  and  of  zero  longitudinal  normal 
force.  In  case  of  stiffening  by  stringers,  the  longitudinal  moment  in  the  boundaries  is  also  not  zero. 
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Hence  in  spite  of  the  requirements  of  Eqs.  (99)  and  (100),  the  displacements  yield 

v  ji  0 

Nx  ^  0  at  x  =  1  and  x  =  x2  (101) 

M,  £  0 

In  order  to  evaluate  the  effect  of  this  non-compliance  with  boundary  conditions  upon  the  critical 
load,  modified  boundary  conditions  were  proposed  in  Ref.  1.  The  ends  of  the  shell  were  assumed  to  be 
elastically  restrained.  For  a  stiffened  conical  shell,  with  such  modified  boundary  conditions,  the  for¬ 
ces  and  the  moments  appearing  in  the  elastic  supports  are  equal  to  the  internal  forces  and  moments 
acting  at  the  boundaries.  The  restraining  forces  are  obtained  from  Eqs.  (43).  At  x  =  1  and  x  =  x2 

Nx  =  [Eh/(1  -v2)]  [(1  +  nl/x)\i'X  +  v  (v^/x  sin  a  +  u/x  -  w  cot  a/x)  -  ^  w(XX/x]  (102) 

=  [Eh/2  (1  +  v)]  [v>x  -  v/x  +  u^/x  sin  a]  (103) 

Mx  =  -lD/a)[(l  +  >701/x)  w>xx  +  v(w>x/x  +  w^/x2  sin2  a)  -  u^/x]  (104) 


The  spring  constants  of  the  elastic  supports  are  defined  as 


ku  ~  Nx/hau  (psi  per  inch) 
hy  =  Nxgi/hav  (psi  per  inch) 

=  6Mx/h8  w  x  (psi  per  inch) 


at  x  =  1  and 


Hence  for  x  =  1,  and  x  =  x. 


JS^cosfr^toxtfyA^  +n/3Arn(l+/il/x)-i.(tBin/sin  a- Au)- Xlx'2n 

1 _ 


E  (1  —  v2)ax  N 


Cn  cos  (n  In  x)Bin 
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2  Cn  cos(n/3  In  x)[Bin  (y-1)  +  n£  Bfn  +  t  Ain/sin  a] 


^  2  (1  +  v)ax 


1  Cn  cos  (n/3  lnx)Bio 

D  =  1 

N 

2  Cn  cos  (n  /3  In  x)[n  /3  (2y-  l)(l+ijw/x)  +  wn  /Q-£x(y  AJn  +  n  /3  A(n)] 


E  2(1—  i/2)ax 


£  Cn  cos(n/3ln  x)n/Q 

n  =  l 


In  order  to  compare  the  magnitude  of  the  spring  constants  for  stiffened  shells  with  those  for  un¬ 
stiffened  ones  the  following  cases  are  considered: 

TABLE  |. 


a  =  30° 


Shell  Type 


a  =  57.59” 


h  =  0. 1  ” 


x2  =  1.5 


\/boh  ej/h  12Iu/boh3  A2/aoh  e2/h  12I22/aoh» 


a.  unstiffened 


b.  stiffened  by  internal  frames 


c.  stiffened  by  external  frames 

d.  stiffened  by  frames  (with  their 
eccentricity  neglected) 


e.  stiffened  by  internal  stringers  0.1471 


f.  stiffened  by  external  stringers  0.1471 


Taking  v  =  0.3  and  n  »  1,  one  obtains  at  i  =  1  (where  the  largest  values  occur): 
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TABLE  II. 


Shell 

a 

b 

c 

d 

e 

f 

VE 

-  0.0059 

-  0.0056 

-0.0071 

-  0.0061 

-  0.0056 

-  0.0048 

Kv/E 

-  1.43 

-  2. 16 

-  1.13 

-  1.57 

-1.42 

-  1.71 

K  /E 

m  ' 

- 

1 

- 

-  0.0052 

+  0.0027 

The  comparison  shows  that  the  magnitudes  of  the  spring  constants  of  unstiffened  and  stiffened 
conical  shells  are  of  the  same  order.  Since  the  values  of  k^  are  very  large,  the  solution  approaches 
the  realistic  boundary  condition  of  zero  circumferential  displacement  (v  =  0).  The  influence  of  ky  upon 
the  critical  load  was  shown  in  Ref.  2  to  be  negligible. 

Another  method  for  evaluation  of  the  effect  of  the  non-compliance  with  the  boundary  conditions  and 
the  resulting  reduction  of  the  critical  load  is  proposed  here. 

In  the  expression  of  the  virtual  work  during  buckling,  Eq.  (61),  the  following  integral  appears. 


27 r 

+  h 


axfN^  aSu 


+  N^^aSv  —  Mx5(w 


„»T* 
'*  X  =  1 


a  d(f> 


(107) 


This  is  the  virtual  work  done  by  the  internal  forces  and  moments  appearing  at  the  boundaries,  or 
the  "boundary  work”.  Since  the  internal  forces  and  moments  at  the  boundaries  are  equal  to  the  forces 
and  moments  acting  upon  the  assumed  boundary  springs,  the  “boundary  work”  is  equal  to  the  work  done 
by  these  fictious  springs.  Now  by  calculation  of  the  critical  load  with  or  without  the  virtual  work  done 
by  the  “springs”,  the  effect  of  every  “spring”  separately,  or  of  all  the  “springs”  together,  can  be 
evaluated.  For  typical  sheila,  stiffened  by  stringers,  the  effect  of  M  (km)  was  found  to  be  less  than 
0.2%.  This  is  not  surprising  since  as  a  result  of  the  definition  of  y,  given  in  Eq.  (56),  the  longitudinal 
moment  at  the  boundaries  is  nearly  zero.  The  effect  of  (ku)  upon  the  critical  pressure  p  is 
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the  most  pronounced,  being  1%  of  pcr  for  short  shells,  and  3%  —  4%  for  long  ones.  It  is  interesting 
to  note  that  the  effect  of  (kv)  is  indeed  small,  as  was  assumed  in  Ref.  1,  being  less  than  1%  for 
both  short  and  long  shells. 

It  has  been  shown  in  Ref.  14  that  by  the  method  of  virtual  displacements  the  equilibrium  boundary 
conditions  are  fulfilled  for  the  complete  displacements,  although  every  term  of  the  infinite  displacement 
series  does  not  fulfil  them.  Hence,  the  magnitude  of  and  N  at  the  boundaries,  for  every  term,  can 
only  affect  the  rate  of  convergence,  and  not  the  solution  itself,  if  their  "boundary  work”  is  taken  into 
account. 

The  method  of  virtual  displacements  requires  that  the  displacements  fulfil  the  geometrical  boundary 
conditions.  The  non-compliance  of  the  condition  of  zero  circumferential  displacement  at  the  supported 
ends  has  therefore  to  be  looked  into.  It  should  be  noted  that  although  v  does  not  vanish  at  the  bounda¬ 
ries,  it  is  nearly  zero  there.  In  order  to  show  this,  the  circumferential  displacement  of  a  typical  shell, 
shell  (b),  has  been  computed.  Minimum  p£f  appears  when  the  number  of  circumferential  waves  t,  is 
10.  Then 

Vj  =  cos (10 <£)  x0-35  [-0.0818744  sin  (ft  lnx)  +  0.00065850  cos  (/3  In  x)]Cj 

v2  =  cos(10<£)x°-35  [-0.07556033  sin  (2/3  In  x)  +  0.00010874  cos  (2/3  In  x)]  C2  (108) 

From  the  above  equations  it  can  be  seen  that  at  the  boundaries  is  only  about  1/100  of 

Vma*  appearing  in  the  shell.  Nevertheless,  the  exact  effect  of  the  displacement  v  at  the  boundaries 
upon  the  critical  load  can  be  estimated  only  by  consideration  of  the  "boundary  work”  of  N  Such  an 
analysis  shows  that  the  effect  is  less  than  1%. 

It  may  be  concluded,  therefore,  that  the  boundary  conditions  applied  here  for  circular  conical 
shells  stiffened  in  two  directions,  differ  only  slightly  from  the  usual  simple  support  conditions,  and 
Eqs.  (99)  and  (100)  are  satisfied  with  sufficient  accuracy. 


13.  COMPLIANCE  WITH  BOUNDARY  CONDITION  OF  ZERO  CIRCUMFERENTIAL  DISPLACEMENT 


The  effect  of  approximate  compliance  only  with  the  boundary  condition  v  =  0  on  the  critical 
pressure  is  very  small.  However,  if  required  (for  example  for  some  other  type  of  load),  it  is  possible  to 
fulfil  this  boundary  condition  for  the  complete  displacement  although  not  for  every  term  of  its  series. 
Then, though 

vn  =};  0  at  x  =  1  and  x  =  x2  (109) 

one  can  prescribe 


v=£  vn  =  0  at  x  =  1  and  x  =  x2 


From  Eqs.  (81)  the  displacement  v  at  the  boundaries  is  : 


vn  =  cos  t  <f>  B.n  cos  (n  /9  In  x)  Co  at  x  =  1  and  x  =  x2 


1  Bin  cos  (n  /3  In  1)  «  B.n 
x2  B.n  cos  (n  0  In  x2)  =  Bjn 


one  obtains 


v  =  cos  t  d>  B.  C  at  x  =  1 

n  ~  in  n 


Vn  =  cos  t  4>  Bjn  cn  at  X  =  x2 


and  the  boundary  condition,  Eq,  (110)  becomes 


S  B.  C  =0 

n  =  I  in  n 


(113) 
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Bj  n  Cn  =  0  (114) 

F.qs.  (114)  represent  the  necessary  connections  between  the  coefficients  C  for  the  compliance 
with  the  boundary  condition  of  zero  circumferential  displacement.  Hence,  now  not  all  the  coefficients 
Cn  are  arbitrary.  It  may  be  assumed  that  the  last  two  coefficients  CN_j  and  CN  are  functions  of 
the  other  N  —  2  arbitary  ones.  The  solution  of  Eqs.  (114)  yields  then 

N—  2 

CN-1  =  g-  Cn  Whefe  gin  =  gin  <Bin  -  Bjn) 


CN  =  n~,  ein  Cn  where  gjn  =  gjn  (Bln  ,  BJn) 


and  the  radial  displacement  can  be  written  as 

N-2  N-2  N-2 

w  =  sin  t  (jj  |  £  x  sin  (n  ft  In  x)  Cn  +  x  sin  [(N  —  1) y3 1  n  x]  £  +x  s*n(N  /3  In  x)£  g.  CB 

n  =  l  n  — 1  n  n  =  l 


Since  one  can  also  obtain  u  and  v  in  a  similar  manner,  the  displacements  become 

N“2  yf 

u  =  sin  t  <f>  £  x  ^  [Arn  sin  (n  ft  In  x)  +  A.n  cos  (n  ft  In  x)] 

+  Sin  '  Ar(N_i)  sin  [(N— l,/3ln  x]  +  Ai(N  l)  cos  [  (N —  1 )  /3  In  x]i 
+  gjn  t  ArN  sin  (N  ft  In  x)  +  A.n  cos  (N  ft  In  x)]j  Cn 

N-2  yf 

v  =  cos  t  0  £  x  1  [  Bfn  sin  (n  ft  In  x)  +  B.n  cos  (n  ft  In  x)] 

+  8in  *  Br(N-i)  sin  WN-1)  ft  In  x]  +  B,(N  i)  cos  [(N-l)£  In  x] ) 

+  gin  sin  (N  ft  in  x)  +  B.n  cos  (N  ft  In  x)]jcn 
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N"2  v 

w  =  sin  t  4>  2^  x  |  sin  (n  $  In  x)  +  gJn  sin  [(N  —  1)£  In  x]  +  g.n  sin  (N  /3  In  x)|  Cn  (116) 

Where  all  the  N— 2  coefficients  Cn  are  arbitrary. 

The  displacements  given  by  Eqs.  (116)  fulfil  all  the  geometrical  boundary  conditions  of  simple 
supports.  Every  term  of  the  displacements  fulfils  the  boundary  condition  of  zero  displacement  in  the 
radial  direction  (w  =  0),  and  although  not  every  term  fulfils  the  boundary  condition  of  zero  circumferen¬ 
tial  displacement  (v  =  0)  the  whole  series  does. 

By  substitution  of  Eqs.  (43)  and  (116)  into  Eq.  (63),  and  integration  (with  the  “boundary  work” 
included)  one  obtains,  since  the  first  N— 2  coefficients  Cn  are  arbitrary, 

N_2 

j*  Cn  t  T(n,m)  +  g.n  T  (N- 1,  m)  +  g^  T (N,  m) 

+  8,  T(n,  N-l)  +  gin  glm  T(N-l.N-l)  +  gjn  giB  T(N.N-l) 

+  8im  T(“.  N)  +  gln  gjm  T(N-  1,  N)  +  gjn  gjm  T  (N.N)]  =  0 

(m  =  1  to  N  -  2)  (117) 

where  T(n,m)  is  the  virtual  work,  multiplied  by  some  known  coefficient,  n  denotes  the  stress  state 
and  m  the  state  of  deformation.  The  lowest  eigenvalue  of  the  determinant  of  the  coefficients  of  CQ 
of  Eqs.  (117)  yields  the  critical  pressure  for  the  case  when  also  the  v  =  0  boundary  condition  is  ful¬ 
filled. 

The  above  solution  can  be  expressed  in  a  more  general  form  with  the  aid  of  Lagrangian  multi¬ 
pliers. 

Since  now  not  all  the  coefficients  Cn  are  arbitrary,  substitution  of  Eqs.  (43),  (81)  and  (82)  into 
Eq.  (63)  yields 

N  N 

2  2  C  T  (n,  m)  8  C  =0 

n=l m=l  n  '  ’  ' 


(118) 
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To  fulfil  the  boundary  condition  v  =  0,  the  coefficients  Cn  and  their  variations  must  be  related 
according  to  Eqs.  (114),  or 


SB.  SC  =0 

m—  I 


S  B.  8  C  -  0 

_  l  J  ru  m 


m  =  l 


(119) 


Multiplication  of  Eqs.  (119)  by  Lagrangian  multipliers  and  addition  to  Eq.  (118)  yields 
N  N  N  N 

S  1  C  T  (n,  m)  5C  +  A,  S  B.  8  C  +  A.  S  B.  8C  =0 

n  — 1  _ ,  n  \  /  ra  l  1  m  m  j  .  j  m  m 

n  —  l  in  —  l  m  =  1  m  =  l 


(120) 


and  then 


(121) 


By  proper  choice  of  the  Lagrangian  multipliers  At  and  A.  the  coefficients  of  SCN  |  and 
8  CN  will  vanish.  Since  the  remaining  variations  5  Cm  are  independent,  their  coefficients  must  also 
vanish.  Thus,  one  may  consider  all  the  variations  8  in  Eq.  (121)  arbitrary.  Hence,  Eqs.  (121) 
and  (114)  yield 


SJ  cn  T  (n,  m)  +  A,  B.m  +  A,  BJm  = 


SB.  C  =  0 

n=l  ln  n 


N 

S  B,  C  =  0 


n=l 


jn  n 


(m  =  1  to  N) 


(122) 


These  are  N  +  2  linear  algebraic  equations  with  n  +  2  unknowns:  Cb  (n  =  1  to  N)  A,  and  Aj . 
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Again,  the  lowest  eigenvalue  of  the  determinant  of  the  coefficients  of  Cn  ,  A,  and  A} ,  of  Eqs.  (122) 
yields  the  critical  pressure,  for  a  shell  which  fulfils  also  the  boundary  condition  of  zero  circumferential 
displacement. 


14.  S  0  L  U  T  I  0  N 

By  substitution  of  Eqs.  (43)  into  Eq.  (63)  and  taking  into  account  that  the  first  two  "corrected* 
stability  equations  are  satisfied  by  the  displacements  Eqs.  (81),  one  obtains 

2ff  *2 

I  /l[Eh/a(l-v2)][MlttiM/x-^1UiMA1-v1Wt„,/*  +  v1xw>M8/kaa+)faW(1/xa 

“ *2vV/k3x  +  X2  w,^/x3  8»n2  a  -  \2  w^/k4  x2  sin2  a]5  u 

+  [Eh/a(l— „2)][_X2  w^/x38in3a  +  x2  w>(|6^/k5  x2  s,n3a-y2  w  x(j6/x2sin  a  +  x2  w >x<6/k6x sin  a]5v 
“(D/a3)fw,„„+2w|Hi/i  -  w^/x2  +  wfi/x3  -2  w^x3  sin2  a  +  2w(ij#/x2  sin2  a 

+  4  w,<^./x48in2  a+v/,<H>cM/**sni  a)  +  cot  a  12(a2/h2)(-vux/x  -u/x2  -vyx2sina 

+  cot  a  w/x2)  +  n2  cot  a  12(a2/h2)  (-u/x2  -  v  ^/x2  sin  a  +  cot  a  w/x2) 

+  *»01  w,a***/x  +  £i  (-U,x1*/X)  +  n02  sin4  a  +  2  w^/x4  sin2  a 

+  w,/x  —  w>xx/x  )  +  ^  (u^/x2  -  u/x3  -  1^0^/x3  sin2  a  +  v^/x2  sin  a 
-  v^/x3  sin  a  -  v  ^/x3  sin3  a  +  cot  a  w/x3  +  2  cot  a  w>(^/x3  sin2  a)  +  ^  (w^^x3  sin2  a 
~  2  w,^/x4ain2a  +  2  w,^,/x5sin2a)  +  r,2t  (w(Ix<^/*»  sin2  a  -  W(X<^/x3sin2  a 
+  w  <Ub’/xi  8‘n2  d  +  iPa3  tan  a/D)  (x  w  /2  +  w 

,rr  txx  t* 

+  w,<^/ x  sin2  a)]  S  w  la3  x  sin  a  dxd  ,<£ 
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2  n 

+  /  |  [Eh/(  1  —  v2)  ]  ax  [(1  +  /ij/x)U|X  +  v(v^/x  sin  a  +  u/x  -  cot  a  w/x)  -  w  xx/x]a5u 
0 

+  [Eh/2(1  +  v)]  ax  [v  x  -  v/x  +  u^/xsina]afiv 

x=*2 

+  Dx[(l  +»701/x)w>xx  +  *(wjX/x  +  v  w^/x2  sin2  a)-  ^  UjX/x](Sw)>x  |  sinad./.  (123) 

Equations  (81)  and  (82)  are  then  substituted  into  Eq.  (123)  .  Since  all  SC  are  arbitrary  one 
can  take  all  SCn  to  be  zero  except  one,  S  Cm.  Then,  after  integration  and  division  by  a  constant, 
rrDsina,  one  obtains 

N 

£  Cn  t  12 (a/h)2  [RC1'(n,  ir.)  +  RC2’(n,  m)  +  RU>  (n,  m)] 

n  =  l 

+  Kl  R>?01(n>  ra)  +  Cl  (n,  m)  +  nlt  R^1'  (n,  ra)  +  C2  R^*  (n,  m) 

+  r?2,  R?,2t  (n,  m)  +  R0T,°2(n,  ra)  +  12  (a/h)2  cot  a  Rh/i2(n,  ra)  +  Ap  RP  (n,  m)] 

+  [12  (a/h)  R  (n,  ra)  +  6  ( 1  —  v)(a/h)  R  ^(n,  m)  +  R  M*(n,  m)]  |  =  0 

(m  =  1  to  N)  (124) 

where  all  the  terms  are  non-dimensional, 

A.p  =  p  a3  tan  a/D 

and  D  =  Eh3/12  (1  -  v2)  (125) 

The  first  part  of  Eq.  (124)  represents  the  virtual  work  done  by  the  “error-loads”  of  the  first  two 
stability  equations.  The  expression  included  in  the  second  square  brackets  represents  the  virtual  work 
of  the  third  stability  equation,  and  the  expression  in  third  Square  brackets  -  the  "boundary  work”.  The 
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N* 

superscripts  of  the  R's  indicate  which  part  of  the  virtual  work  they  represent.  For  example,  R  (n,m) 
represents  the  virtual  work  done  by  the  longitudinal  force  N  at  the  boundaries. 

If  in  Eq.  (124)  the  expression  included  in  the  curled  brackets  is  denoted  T(n,m),  the  equation 
may  be  written  as 

Cj  T  (1,1)  +  C2  T  (2,1)  + . +  CN_i  T  (N-1,1)  +  CN  T  (N,l)  =  0 

Cj  T(l,2)  +  C2  T  (2,2)  + . +  CNl  T  (N-1,2)  +  CN  T(N,2)  -  0 


Cj  T(1,N-1)+  C2T(2,N-1)  + 
CjTO.N)  +  C2T(2,N)  + 


+  CN_1  T(N-1,N-1)+  CnT(N,N-1)  =  0 
+  cn_it<n-1-n>  +  CN  T  (M,  N)  =  0  (126) 


The  lowest  eigenvalue  of  the  determinant  of  the  coefficients  of  Cn  yields  again  the  critical 
pressure  for  general  instability 

|  T  (n,  m)  |  =  0  (127) 

The  integral  value  of  t  (the  number  of  circumferential  waves)  which  makes  pcr  a  minimum  must 
be  used  in  calculations. 

The  R  functions  of  Eqs.  (124)  are  defined  as  follows  : 


(1)  2V  —  1  1A  IB  2v— l  1A  IB 

R  (n,m)  -  1,  (n,m)  [Arm  Fj  (n)  +  Brm  Fj  (n)]  +  i/  (n,  m)  [Arm  F2  (n)  +  Brm  Fa  (n)] 


+  I2  (m,n)[Aiin  F|A(n)  +  pj®  (n)]  +  \n.rn)  [A,m  F12A(n)  +  Bim  F2 
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R<2)(n,m)  =  I ly  2(n,m)[Arm  fJA(o)  +  Bfm  fJ%)]  +  l*/  2(n,m)[Atm  F2%)  +  Bfm  F2B(n)] 

t  2V— 2  2  A  2  B  » y  o  ai  2H 

+  l2  (m,  n)  tAlm  Fj  (n)  +  B.m  Fj  (n)]  +  l/  (n,m)[Ai-F2  (n)  +  BimF2  (n)] 
R<3)(n,m)  =  ljy_3(n,m)  Afm  FjA(n)  +  I2y_3(n,m)  Afm  fJA(o) 

+  3(m.n)  Alm  FiA(«)  +  Aim  F3A(n)  O' 


R*70  (n,m)  «  I2y  4(n,m)  Fj01(n)  +  I2y  4(n,m)  F201(n) 
R^1  (n,m)  -  I2y  3(n,m)  F'f1  (n)  +  I2y  3  (n,m)  F^(n) 
R,,,1(nIm)=  I2y  4(n,tn)  FjU(n)  +  I,2y  4(n,m)  F^n) 
R^2(n,m)  =  I2y  2(n,tn)  F^(n)  +  I2y  2(n,m)F22(n) 
r7?  (ni  m)  =  3(n,m)  Fjt2  (n)  +  I2y  3(n,  m)  F2‘2(n) 


0l?02  2y-i  .  ,,01)02  2y—3  0)?02 

R  (n,m)  =  Ij  (n,  m)  Fj  (n)  +  I2  (nfm)F2  (n) 


RhM2(n,m)  =  I2y  \n,  m)  F^%)  +  I2y  ^n,  m)  F^2  (n) 


RP(n,m)  =  I3y(n,m)  FP(,n)  +  I2y  (n, m)  FP  (n) 


_Nx  <vm+n  2 y  _Nj  .  „N» 

R  (n,m)  =  Alm  [(-1)  x/F2  (n)  -  Fj  (n)] 


4  .m+n  2y~  2  Mx  _Mx 

R  (n,m)  =  m/3  [(-D  x2  F2  (n)  -  F.A  (n)] 


(130) 
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FilA(n)  -  Acn  -  v12  a3n  +  *23y  -  *24  t2/ain2  a 
Fi1B(n)  =  X2S  t3  /sin3  a  -  x26  yt/sin  a 

F2  (D)  =  ^11  Adn  —  Xl2  b3n  +  X23  11  $ 

F2  (n)  =  “  *26  11  ^  t//8‘n  a  0^1) 

2  A 

F!  («0  -  -  /*!  Acn  -  X2  Y  +  X2  t2/sin2  a 

2B 

F!  (n)  =  -  x2  t3  /  sln3  a  +  X2  y  t/sin  o. 

F 2  (°)  =  “  /*!  Adn  -  X2  “Z3 

Fa  (n)  =  +  \2  n  /3  t/sin  a  (132) 

F1  (“)  =  Xx  a3n 

F2A(“)  =  X,  b3n  (133) 

FJ01(n)  =  y  (y— 1)  (y-2)  (y-3)  +  n2  /S*  [n2  p2  -  11  +  6y  (3-y)] 

Fj01  (n)  =  2n/3  (n2  jQ2(3-2y)  +  2  y3  -  9y2  +  11  y  -  3]  (134) 

Ff(n)  =  -  A  n  a3n  +  AJn  b3n 

F  2 1  (n)  “  -  Arn  b3  □  "  Ain  a8n 


(135) 
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F^  (n)  =  (t2/sin2  a)  (-y2  +  3y  —  2  +  n2/9) 

F^n)  =  (t2/sin2  a)  (3  -  2y)a/9  (136) 

A 

Fj  (n)  «=  Afo  ( y  +  t2 /sin2  a- 1)  +  Ain(-n|S)  +  Btn  [(t/sin  a)  (1-y)  -  t8/sin3  a] 

+  BJn  nj9  (t/sin  a)  +  cot  a  (1  -  2t2/sin2  a) 

& 

F2  (n)  =  Afn  n/3  +  AJn  (t2/sin2  a  +  y  -  D  +  Bfn  n  /9  (-t/sin  a) 

+  Bln  [- t3/sin3  a  +  (t/sin  a)  (1  -  y)]  (137) 

(n)  =  (t2/sin2  a)  (- y2  +  2y  —1  +  n2  fi2) 

F2  (n)  *=  (t2/sin2a)  (2  -  2y)n0  (138) 

F^  (n)  ■=  (t4/sin4  a)  (1  +  ?j02)  +  (2t2/sin2  a)  (n2 /92  -  rj02  -  2  +  2y  -  y3)  +  n4/9* 

+  “2  P  (>?0 2  ~  4  +  12 y  -  6y2)  +  *j02  (2  y  —  y2)  +  4y2  -  4y3  +  y4 

0770  2  2 

F2  (n)  =  4n$  [(t2/sin2a)  (1-y)  +  a2  /3  (1-y)  +(rj02/2)  (l~y)+ 2y  -  3y2  +  y3]  (139) 

hfX2 

Fj  (n)  “  Arn  C  — *^y  ~  O  +  ^  +  Ai„l'n0  +  B,n  (t/sin  a)  (1  +  H2)  +  COta(l  +  fi2) 

F^2(n)  =  Arn  (-  ynfi)  +  A.q  [~vy  -  (1  +  n2)]  +  Bin  (t/sin  a)(l  +  /j2)  (140) 

Fj  (n)  =-  t^/sin2  a  -  (1/2)  (n2  /Q2  —  y  —  y2) 

F^(n)  -  (l/2)n/3(l  +2y)  (141) 
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F^X(n)  «*  (A lDy  +  Afn  n /9)  (1  + /ij) -v[B.n  (t/sino)  -  Ajn]  +  n^(l  -  2y) 

=  (Aln  y  +  \D  n$)  (1  +  gj/x2)  -v[Bin  (t/sin  a)  -  A.J  +  np  (1-2 y)/x2 

F  X^(n)  =  -  Bln  (1  -  y)  +  Brn  n  /3  +  A,n  (t/sin  a) 

F^X(n)  =  nfi  [(2 y  -  1)  (1  +  i)Q1)  +  v]  -  £,  (A.n  y  +  \B  n  0) 

Mr 

F2  (n)  =  n/9  [(2y- 1)  (1  +  ?j01/x2)  +  v]  -  ^  (A in'/+  Arn  n /3)  (142) 

15.  APPROXIMATE  SOLUTION  BY  NEGLECT  OF  ECCENTRICITY  OF  THE  STIFFENERS 

For  conical  shells,  stiffened  by  equal  and  equally  spaced  frames,  a  simple  approximate  method 
for  calculation  of  the  critical  pressure  can  be  derived.  Calculations  of  the  critical  pressure  for  ring-stif¬ 
fened  conical  and  cylindrical  shells  (see  Ref.  13),  have  shown,  that  when  the  eccentricity  of  the  stif¬ 
feners  is  neglected,  pcr  is  between  that  for  internal  rings  and  that  for  external  ones.  Hence,  for  the 
purpose  of  an  approximation,  the  following  assumptions  may  be  made  for  a  ring  stiffened  shell :  the 
extentional  stiffness  of  the  shell  is  increased  by  that  of  the  rings,  and  the  total  bending  stiffness  in 
the  circumferential  direction  is  taken  to  be  the  stiffness  of  the  combined  cross-section.  The  torsional 
stiffness  of  the  ring  is  assumed  to  be  small  and  hence  negligible.  With  these  assumptions,  the  inter¬ 
nal  forces  and  moments  given  in  Eqs.  (36)  and  (37)  become  : 

Nx  =  [Eh/(1  -  v2)][<x  +  v^] 

=  [Eh/(1  -  v2)][(l  +  (r2)  ^  +  v<x] 

,0  =  =  [Eh/2  0  +  v)]  y^ 


N 


(143) 
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Mx  =  ~  (D/a)  (kx  +  vk^) 

U4>  =  -(DAHd  +  S2t?02)V  + 


-(D/a)(l-v)^ 


(144) 


If  one  denotes 


^2  rlo2  ^2 


(145) 


and  equates  the  circumferential  bending  stiffness  given  in  F.qs.  (144)  with  that  of  the  combined  cross- 
section,  one  obtains 


D  ( 1  +  >?2)  -  (l/a0)[E2I22+  E,  A2(e2  -  Z2)2  +  a0D  +  EaQh  z'/(l-v2)] 


(146) 


and  hence 

r)2  =  12(1  -  v2)  (E2/E)  II22/a0h3  +  (A2/ah)  f(e2  -  E2)/h]2  I  +  12  (E2/h)Z  (147) 

With  the  simplified  force  and  moment  expressions  of  F.qs.  (143)  and  (144)  no  “correcting  coeffi¬ 
cients”  are  needed.  In  all  the  formulae,  the  terms  introduced  by  the  stiffeners,  except  /j2  and  rjQ2  , 
vanish  ;  and  of  the  two  exceptions  r)Q2  is  replaced  by  tj2  from  F.q.  (147). 

A  further  simplification  is  possible  if  one  neglects  ^2.  Bodner  (Ref.  9)  showed  on  mathematical 
grounds  that  neglecting  of  n2,  introduces  an  error  of  less  than  1  percent  in  pcr.  Calculations  for  ty¬ 
pical  cylindrical  and  conical  shells  carried  out  in  connection  with  the  present  work  verified  that  the 
error  is  much  smaller  than  1  present  in  both  types  of  shells.  Furthermore,  the  calculations  showed, 
that  the  approximation  involved  by  neglecting  ^  is  of  much  smaller  magnitude  than  the  neglect  of  the 
effect  of  eccentricity  of  stiffeners  inherent  in  the  approach  of  this  section. 

Hence,  in  the  approximate  solution  proposed  here  for  ring-stiffened  conical  shells  under  hydro¬ 
static  pressure,  the  following  substitutions  must  be  made  in  all  formulae  : 


so 


f2! 

-  0 

f*2 

0 

*1 

-  0 

*2 

-  0 

*1 

->  0 

4 

0 

-  0 

*01 

->  0 

*.2 

0 

*02  "*  *2 


(148) 


With  Eqs.  (148)  the  calculation  of  pcr  becomes  much  easier,  since  all  the  terms  involving 
“correcting  coefficients”  are  multiplied  by  /r^,  or  which  are  assumed  to  be  zero. 

The  final  results  of  the  above  approximate  solution  were  also  obtained  (at  the  same  time),  inde¬ 
pendently,  in  Ref.  10,  by  another  approach.  The  method  applied  there,  was  substitution  of  the  stiffened 
conical  shell  by  an  equivalent  orthotropic  one.  One  very  small  difference  should,  however,  be  noted  . 
in  Ref.  10  the  increase  of  the  shell  cross-section  due  to  the  rings  (when  it  is  taken  into  account)  is  not 
multiplied  by  (1  -  v2)  whereas  here  /i2,  which  stems  directly  from  Eqs.  (38), contains  this  multiplica¬ 
tion  already. 

Although  the  same  final  results  are  obtained  by  the  equivalent  orthotropic  shell  approach  and  by 
that  of  this  section,  the  two  approaches  differ  considerably.  The  orthotropic  approach  requires  ortho¬ 
gonal  shell  properties,  which  imply  equal  and  equally  spaced  rings  and  stringers  (if  any)  which  vary 
according  to  the  cone  radius;  whereas  in  the  approximate  method  of  this  section,  as  well  as  in  the 
preceeding  more  accurate  method,  no  similar  restrictions  are  implied. 

Here  again,  the  additional  stiffnesses  due  to  the  stiffeners  can  be  some  functions  of  x  (the 
stiffeners  have  to  be  symmetrical  with  respect  to  the  cone  axis). 

If  /r  and  pt2  and  the  eccentricity  of  the  stiffeners  are  neglected,  the  internal  forces  and 
moments  become 


Nx  =  [Eh/(1  -v2)](<x  +  v<^) 
«  [Eh/(1  -  v2)]  (^  +  w.) 
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Nx^  =  V  =  [Eh/2(1  +  V)]}/^  <149> 

Mx  =  (-D/a )|kx  [1  +  r)l  (x)]  +  VK^  ! 

U(j>  =  (-D/a^K^Ll  +  V2  (x)l  +  VKx  * 

M*^6  =  ~  M«^x  =  (D/a)  0  ~  ^  Kx<f>  (ISO) 

It  should  be  noted  that  rjl  is  a  function  of  x  even  when  the  cross  section  of  the  stringers  does 
not  vary.  For  constant  area  stringers,  one  obtains 

D  [1  +  7?!  (x)]  =  (l/b0x)  [Fj  Iu  +  Et  At  (ej  -  Zj)  +EbQ  xhz^l  -  v2)  +  bQ  xD]  (151) 
and  then 

ijl(x)  =  (l/bQ  Dx)  (F1  Iu  +  Ft  Aj  (et  _  zt)2  +  E  bQ  x  h  z\  /'(l  -  u2)  ]  (152) 

From  Eqs.  (36)  one  obtains  directly 

»?i(x)  -=  7?01/x  =  T)l/x  (153) 

Two  similar,  but  not  equal,  functions  for  rj j  are  obtained.  The  expression  in  the  square  brac¬ 
kets  of  Eq.  (152)  varies  only  slightly  with  x.  For  example,  for  A;/b0h  =  0.1471,  e^/h  =  1.653, 
Ijj/boh3  =  0.7819  and  E1  =  E  this  value  is  7%  larger  than  that  at  x  =  1  when  x  =  6,  or  11% 
larger  when  x  oo .  rjx  (x)  defined  in  F.q.  (153)  is  the  more  accurate  function,  since  it  stems  directly 
from  Eqs.  (36).  Eqs.  (152)  and  Eq.  (153)  can  be  exactly  equal  only  at  a  particular  value  of  x,  say  at 
the  midheight  of  the  conical  shell,  where 

x  =  (1  +  x2)/2  (154) 


then  rj  1  becomes 
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Vl  =  12  (1  —  v2)  (Ej/E)  i  Ij  1/b()h3  +  (Aj/b.jhjKe,  _  5  t)/h]2  1  +  12^ /h)2  i  (155) 

For  a  conical  shell  stiffened  by  stringers  only,  all  the  terms  due  to  the  stiffeners  vanish,  except 
1701  which  must  be  replaced  by  jjj  from  Eq.  (155). 

Since  the  effect  of  stringers  on  p  is  rather  small,  and  the  effect  of  their  eccentricity  is  of  the 
same  magnitude,  the  calculation  of  pcr  without  consideration  of  the  eccentricity  has  little  value. 

Hence  Eq.  (155)  will  be  of  limited  use  for  stringer-stiffened  conical  shells  under  hydrostatic  pressure. 
It  may  however  be  usefull  in  other  loading  cases. 


16.  EFFECTIVE  LENGTH  OF  SHEET 

In  the  previous  sections  it  has  been  assumed  that  the  stiffeners  are  closely  spaced  and  therefore 
the  entire  shell  is  active.  If  the  distance  between  the  stiffeners  is  larger  than  a  certain  magnitude, 
only  part  of  the  sheet  between  the  stiffeners  is  active.  The  resulting  decrease  in  the  total  stiffness 
of  the  shell  usually  expressed  as  “effective  length  of  the  sheet”  can  also  be  expressed  as  a  decrease 
in  the  modulus  of  elasticity  of  the  stiffeners. 

If  the  circumferential  stiffness  of  the  combined  cross-section  for  a  wholly  active  sheet  but  “ef¬ 
fective  moduli"  of  stiffeners  is  compared  with  that  for  an  “effective  length”  —  ae  -  one  obtains . 

E2  [I22  +  A2  (e2  -  +  E  Kb3/12*1  -  y2)  +  a0h  h  A1  -  u2)]  =  E2  [I22  +  ^2  “  *2^ 1  + 


+  E[a  h3  /  12(1  —v2)  +  a  h  z';2/(l  -  v2)] 


•  2 


(156) 


and  then 


e2/e  = 


12(1— v2)  (E;/E)|I22/aoh3  +  (A2/aQh)[(e2  ~z^hf\  +  ae/aQ-  1  +  12(a»2/h)%e/a#)-  12 (z 2/h)2 


(157) 


12(1  -  v2)  1 122/aQh3  +(  A2  /aQh)  [(e2  -  ijh  ]'  } 
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In  the  longitudinal  direction  it  is  assumed  that  the  "effective  length”  varies  lineary  in  the  same 
manner  as  the  actual  distance  between  the  stringers.  The  respective  cross-sections  are  equated  at 
the  midheight  of  the  conical  shell. 

Then 

12  (1  _^2)(E\/E)U11/b0h3  +  (Aj/b^ttej-SJ/h  fl  +  S[b  A>0-  1  +  12(5»1/h)2(be/b0>-12(S1/h)2] 


VE  = 


12(1  — v2){I11/b0h3  +(A1/b0h)[(e1-z])/h  f\ 


(158) 


In  Eqs.  (156)  to  (158),  ao  and  be  are  “effective  lengths”,  Ej  and  E2  are  “effective” 
moduli  of  elasticity  of  the  stiffeners,  E*j  and  E*2  are  actual  moduli  of  elasticity  of  stiffeners, 
and  t2  are  the  distances  of  the  overall  centroid  of  the  stiffener-shell  combination  from  the  middle 
surface  when  the  sheet  length  is  or  b^x  respectively,  and  and  z’  are  the  distances  of  the 
overall  centroid  of  the  stiffener-shell  combination  from  the  middle  surface,  when  the  sheet  length  is 
ae  or  bex  respectively. 


17.  APPROXIMATE  FORMULAE  FOR  RING  STIFFENED  SHELLS. 

In  Ref.  10  it  is  shown  that  the  ratio  between  the  critical  pressure  pcr  of  an  orthotropic  conical 
shell  and  that  of  its  equivalent  cylindrical  shell,  pcr  may  be  approximated  by  the  same  ratio  for  iso¬ 
tropic  shells. 

Pc /Per  = 

=  1  -  (Rt/R2)  =  1  -  (l/x2)  (1S9> 

The  function  g  (</i)  is  given  in  Ref.  1  and  Ref.  11.  It  is  reproduced  here  in  a  tabular  form  (Table  3), 
the  values  having  been  read  off  the  curve  in  Ref.  11,  It  was  pointed  out  in  Ref.  10  that  the  accuracy 
of  g  (i p)  diminishes  slightly  when  a  is  greater  than  45%. 
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TABLE  3 


0 

0 

0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

g 

1.00 

1.005 

1.01 

1.02 

1.04 

1.075 

1.13 

1.19 

1.22 

1.21 

1.175 

The  equivalent  cylindrical  shell  is  taken  as  one  having  a  length  equal  to  the  slant  length  of  the 
cone,  £  ,  a  radius  equal  to  its  average  radius  of  curvature  pav,  and  the  same  thickness  and  ring- 
stiffeners.  Based  on  the  results  of  Niordson  (Ref.  12)  and  Bodner  (Ref.  9),  the  critical  pressure  for  an 
equivalent  cylindrical  shell  may  be  written 


Pcr/E  =  [l/(tQ  +  0.5  c^)]  Kh/pav)[cX  +  ^  ^  ft3/ 12(1  - 


pL  ][(*; + s)2 + 


0  1 


(160) 


where 


tQ  =  t/cos  a 

pay  »  (a tan  a/ 2)  ( 1  +  %2)  =  (Rt  +  R2)/2  cos  a 
C0  —  rrpa  J-t  =  ( n  tan  a/ 2)[(xa  +  lV£2  -  t)l 

£  =  a(x2  -  1)  =  (R2  -  R  J/sin  a  (161) 


In  general  cQ  is  much  smaller  than 
form  by  changing  the  first  denominator  from 

Pe[/E  -  [l/(to  +  cj)l  !(h/pav)  tc02/(to  + 


tQ  and  then  F.q.  (160)  may  be  written  in  a  more  convenient 

2  2  2  2 
(tQ  +  0.5  c0)  to  (t0  +  cQ) 

+  [h3/12(l-v2)pa3v][(t2+  c2)2  +  T,2tJ)l  (162) 


This  change  reduces  the  critical  pressure  slightly.  Its  physical  interpretation  is  the  conservative 
assumption  that  twice  the  actual  hydrostatic  pressure  acts  upon  the  rigid  bulkheads  closing  the  equivalent 


55 


cylindrical  shell.  Since,  however,  in  case  of  hydrostatic  pressure,  the  effect  of  the  pressure  acting  in 
the  longitudinal  direction  is  small  compared  to  that  acting  radially,  Eq.  (162)  may  be  expected  to  yield 
a  good  approximation. 

One  can  calculate  the  minimum  value  of  pcr  from  Eq.  (162),  by  assuming  p  to  be  a  continuous 
2  2 

function  of  (tQ  +  cQ).  Only  very  slight  conservative  errors  are  involved  in  this  assumption.  Hence, 
with  the  notation 

t02  +  co  =  S  (163) 

One  obtains  for  a  ring-stiffened  conical  shell 

pcr/E  =  [hVp!v12(l-v2)][U  +  n2)S-2n2  c'  +r,2c^/S  +  (pay/h)2  12(1  -v2)  (c0Vs3)]  g  (./,)  (164) 

p  is  a  minimum  when 
rcr 


2 ,0.5  ,0.5 , 


S  “  lo  +  co  =  co  U2  co  +  U2  c0  +  144(! -^23d  +  n2)(pAV/^)  1  I  [2(1+  rj2)l 


Since  cQ  is  usually  much  smaller  than  tQ,  one  could  alternatively  approximate  Eq.  (160)  by 

.  2  ,  .  2 
neglecting  cQ  in  comparison  with  tQ 

Pc/E  ~  OA0)  th  co Aav  to  +  h3  to(1  +  ^2 12 ( 1  -  (166 

If  pcr  is  again  assumed  to  be  a  continuous  function  of  tQ,  it  is  a  minimum  when 


:  =  c0  V6(Pa,A)0-5(l-2)0-25(U,2r°l25  =  (1-^)U'“  ,V6'(p,/)(p.v/hra(l+  ,2)  (167) 


and  this  minimum  value  is 


Pcr/E  =  [  V6n79(l  -  v2)  5  1  (pav/^)  (h/pBV)  (1  +  t)2)  g  (</<) 


(168) 


It  should  be  noted  that  whereas  Eq.  (164)  yields  a  low  value  for  pcr,  a  high  value  is  obtained 
from  Eq.  (168),  since  there  the  effect  of  the  frames  is  slightly  exaggerated. 

4 

Further  simplification  of  Eq.  (164)  is  possible  by  neglect  of  the  small  term,  n2  C0/S-  Th*s 

raises  pcr  slightly.  The  increase  is  however  partly  compensated  by  a  replacement  in  Eq.  (164)  of 

2 

S,  from  Eq.  (16S)  by  its  approximation  tQ  from  Eq.  (167).  The  simplified  formula  obtained  in  this 
manner  and  that  of  Eq.  (168)  ate  averaged,  and  finally  a  simple  formula  for  the  critical  external  hydro¬ 
static  pressure  of  a  ring-stiffened  conical  shell  (equal  rings,  equally  spaced)  failing  by  general  insta¬ 
bility  is  obtained  : 

pcr/E  =  [\f6n/9  (l-i/2)0' 75]  (pav/f. )  (h/ pav)2  5  Kl  +  rj2)°'75-[\/6't/8(l-v2)0'25](pav/£)(h/pBv)a5i ^IgO/O 

(169) 

If  r}2  =  0  is  substituted  in  Eq.  (169),  Seide’s  approximate  formula  for  unstiffened  conical 
shells  (Ref.  11,  Eq.  23)  is  obtained. 

With  v  =  0.3,  sf<on/%  (1  -0-3  )  '  =  0.99~1,  and  Eq.  (169)  becomes 

Pcr/E  =  0.92(Pav/m/Pav)2'5  [U+t?/'75 -(pav/£)(h/pBV)°'%2]  g  W  (17°) 


If  pQ  is  the  critical  pressure  for  an  unstiffened  conical  shell,  which  was  calculated  with 


tj2  =  0,  Eqs.  (169)  and  (170)  may  be  rewritten  in  a  simpler  form,  where  \[(>  -  v2)  is  assumed 

to  be  approximately  unity  for  all  likely  values  of  v , 


Per  =  Po  +  J?2^'75  2 ] 


(171) 


In  this  section,  as  elsewhere  in  this  report,  the  work  of  the  shear  forces  0  is  not  taken  into 
account.  The  theory  is  a  Donnell  type  theory,  and  is  valid  only  when  t  (  the  number  of  circumferen¬ 
tial  waves)  is  greater  than  two  (Ref.  9).  Hence,  the  approximate  formulae  of  this  section  are  applialie 
when  the  i,  obtained  from  Eqs.  (165)  or  (167),  is  more  than  two.  In  the  above  approximate  derivation 
it  was  also  assumed  that  the  semi-empirical  function  for  t  of  an  unstiffened  conical  shell,  when 
i/i  >  0.o4  (Fig.  8  of  Ref.  11)  holds  also  for  a  ring-stiffened  shell. 
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It  may  be  pointed  out  that  all  the  above  approximate  formulae  apply  also  to  the  limiting  case  of 
a  ring-stiffened  cylindrical  shell.  For  the  cylindrical  shell,  pay  of  the  above  formulae  is  the  radius 
of  the  shell,  £  its  length,  and  cos  a  and  g  (tp)  are  unity. 


18.  NUMERICAL  RESULTS  AND  DISCUSSION 

The  critical  pressures  for  general  instability  are  computed  for  typical  cases  (Table  4).  The 
torsional  resistance  of  the  stiffeners  is  neglected  and  the  whole  distance  between  the  stiffeners  is 
taken  as  the  “effective  length"  of  the  sheet.  The  moduli  of  elasticity  of  sheet  and  stiffeners  are  equal. 

Two  cases  are  computed  by  all  the  methods  proposed  in  this  report.  The  results  obtained  are 
compared  with  those  for  unstiffened  conical  shells  given  in  Ref.  1.  Two  other  cases  are  computed 
only  by  the  methods  of  Section  15  (j }Q2  -*  t)2)  and  Section  17  (Eq.  169),  and  compared  with  the  results 
obtained  by  Seide’s  approximate  formula  (Ref.  11)  for  unstiffened  conical  shells. 

The  rate  of  convergence  of  the  solution  is  similar  to  that  for  unstiffened  conical  shells.  For 
short  shells  a  two  term,  or  even  a  one  term  solution  is  sufficient.  For  long  shells,  additional  terms 
must  be  considered. 

The  “boundary  work”  reduces  the  critical  pressure  slightly.  Its  influence  upon  the  critical 
pressure  is  less  than  1%  for  short  shells  (shells  Nos.  2  to  7)  but  it  is  3%  to  4%  for  long  ones  (shells 
Nos.  9  to  12). 

Table  4  demonstrates  clearly  that,  as  for  stiffened  cylindrical  shells,  frames  are  very  effective 
in  stiffening  conical  shells  against  hydrostatic  pressure.  For  short  shells  (shells  Nos.  2  to  5)  ad¬ 
dition  of  only  15%  of  material  increases  the  critical  pressure  more  than  3  times  that  for  the  similar 
unstiffened  shell.  For  long  shells  (shells  Nos.  9  to  12),  this  increase  is  larger,  and  for  the  same 
addition  of  15%  of  material  it  is  more  than  3.5  times.  Addition  of  the  same  material  uniformly  to  the 
thickness  of  the  shell  would  increase  the  critical  pressure  only  1.4  times.  Addition  of  22. 5%  of  material 
(shells  Nos.  14  and  15)  increases  the  critical  pressure  more  than  6  times,  whereas  addition  of  the  same 
material  to  the  thickness  would  increase  the  critical  pressure  only  1.66  times.  The  increase  with  length 
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in  the  ratio  of  the  critical  pressure  for  ring-stiffened  conical  shells  to  that  for  corresponding  unstif¬ 
fened  ones,  indicated  also  by  the  approximate  formula,  is  due  to  a  larger  part  of  the  hydrostatic  pressure 
being  transmitted  in  the  circumferential  direction  for  longer  shells. 

Stringers  are  much  less  effective  as  stiffeners  against  hydrostatic  pressure. 

The  number  of  circumferential  waves,  t,  for  which  a  minimum  of  the  critical  pressure  is  obtained, 
decreases  with  increase  in  the  stiffening  of  the  shell.  The  order  of  magnitude  of  this  decrease  given 
by  Eq.  (167)  and  verified  by  the  results  of  Table  4,  is  ( 1  +  rj2)~°' 125  It  should  be  noted,  that  when 
ip  >  0.64,  the  number  of  waves  obtained  from  Eq.  (167)  is  multiplied  by  a  coefficient  taken  from  Figure 
8  of  Ref.  1 1. 

The  effect  of  the  eccentricity  of  the  frames  on  the  critical  pressure  may  be  summarized  as  follows. 
Internal  frames  yield  higher  general  instability  pressures  than  external  frames.  In  a  typical  case,  inter¬ 
nal  frames  (shell  No.  2)  yield  a  critical  pressure  7%  greater  than  that  obtained  by  external  frames  (shell 
No.  3).  For  long  shells,  the  effect  of  the  eccentricity  of  the  frames  is  much  more  pronounced.  For  a 
typical  long  shell  with  internal  frames  (shell  No.  9)  the  critical  pressure  is  12%  higher  than  for  the  same 
shell  with  external  frames  (shell  No.  10).  This  effect  should  be  taken  into  account  especially  in  the 
analysis  of  experimental  results.  Internal  frames  yield  a  higher  critical  pressure  due  to  their  smaller 
radius  which  makes  them  stiffer.  The  critical  pressure  obtained  by  the  method  of  Section  15  (r702  -»rj2) 
is  found  to  be  somewhere  between  the  critical  pressures  for  internal  and  external  frames. 

For  longitudinal  stiffeners  (stringers),  the  effect  of  eccentricity  is  opposite  to  that  in  frames. 
External  stringers  yield  higher  critical  pressures  than  internal  ones.  In  Table  4,  external  stringers  of 
the  same  magnitude  as  that  of  the  frames  (shell  No.  7)  yield  a  critical  pressure  only  12%  higher  than 
that  for  the  corresponding  unstiffened  shell,  whereas  for  internal  stringers  (shell  No.  8)  the  increase 
is  only  2%.  Note  that  the  inc.ease  in  critical  pressure  due  to  stringers  is  considerably  less  than  that 
obtained  by  uniform  thickening  of  the  shell  with  the  same  amount  of  material.  Hence,  stringers,  and 
especially  internal  stringers  are  very  inefficient  stiffeners  for  conical  and  cylindrical  shells  under 
hydrostatic  pressure  (see  also  Ref.  13).  However,  if  stringers  are  taken  into  account  one  should  not 
ignore  their  eccentricity.  The  above  shown  inefficiency  of  stringers  applies  only  to  the  case  of  external 
pressure  loading.  For  other  loads,  stringers  are  much  more  effective  and  should  be  subject  to  further 
investigation. 
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